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CONVERSION  FACTORS,  NON-SI  TO  SI  (METRIC) 

UNITS  OF  MEASUREMENT 

Non- SI  units  of  measurement  used  In  this  report  can  be  converted  to  SI 


(metric)  units  as  follows: 

Multirlv 

Bv 

To  Obtain 

degrees 

0.01745 

radians 

cubic  feet 

0.02831685 

cubic  metres 

inches 

25.4 

millimetres 

feet 

0.3048 

metres 

feet  per  second  (fps) 

0.3048 

metres  per  second 

inches  per  second  (Ips) 

3.6576 

metres  per  second 

pounds  (mass) 

0.45359237 

kilograms 

pounds  (mass)  per 
cubic  foot  (pcf) 

16.0179 

kilograms  per 
cubic  metre 

pounds  (force)  per 
square  inch  (psi) 

0.006894757 

megapascals 

pounds  (force)  sq'iare  seconds 
per  inch 

175.1265 

kilograms 

pounds  (force)  per  inch 

175.1265 

newtons  per  metre 

EVALUATION  OF  IN-STRUCTURE  SHOCK  PREDICTIQS 
TECHNIQUES  FOR  BURIED  STRUCTURES 

CHAPTER  1 
INTRODUCTION 


1.1  Background 

The  analysis  of  burled  structures  designed  to  resist  blast  loading  Is  a 
continuing  problem  In  the  area  of  engineering  mechanics.  The  complex  struc* 
rural  dynamics  of  the  problem,  coupled  with  the  uncertainties  of  soil  struc¬ 
ture  Interaction,  make  this  a  most  difficult  challenge.  Even  If  Che  structure 
Itself  is  designed  to  survive  Che  expected  dynamic  loads,  Che  shock  environ¬ 
ment  Inside  may  be  severe  enough  to  harm  personnel  or  equipment.  A  thorough 
understanding  of  Che  expected  shock  environment  Is  therefore  needed  to  design 
ways  to  reduce  this  damage.  The  present  simplified  methods  for  the  analysis 
of  this  type  of  problem  are  only  directly  applicable  Co  very  simple  structures 
and  are  thought  to  yield  very  conservative  results  for  more  complex  stxruc- 
Cures.  As  a  result,  large  amounts  of  money  are  now  being  spent  to  mitigate 
shock  levels  which  may  be  unrealistically  high.  This  program  of  study  was 
prompted  by  the  hope  Chat  a  relatively  simple,  fast,  and  accurate  analysis 
procedure  could  be  found,  which  will  give  realistic  estimations  of  In-struc- 
Cure  shock  for  complex  burled  reinforced  cdncrete  structures. 

1.2  Objective 

The  objective  of  this  study  was  to  conduct  and  evaluate  in- structure 
shock  calculations  for  burled  reinforced  concrete  structures  under  d)mamlc 
loads  from  burled  conventional  explosives.  Of  primary  importance  was  the 
evalua*-  of  two-dimensional  (2-0)  finite -element  analysis  techniques  applied 
Co  this  K  three-dimensional  problem.  Two  personal  computer  (PC)  based, 

2-D  flnlte-eiement  programs,  were  evaluated  for  application  to  this  problem. 

1.3  Scope 

This  study  included  simplified  seml-emplrical,  slngle-degree-of-freedom 
(SDOF),  and  2-D  finite -element,  In-structure  shock  analyses  bf  burled  rein- 


forced  concrete  structures,  A  literature  search  was  conducted  to  obtain 
available  data  to  assist  in  the  evaluation  of  these  calculations. 

The  project  began  with  an  analysis  of  the  existing  data  base  of  in* 
structure  shock  tests.  A  recent  series  of  four  dynaoic  tests  conducted  by  the 
US  Aray  Engineer  Vaterways  Experloent  Station  (VES) ,  Structures  Laboratory. 
Structural  Mechanics  Division,  were  selected  as  having  the  best  data  cvirrently 
available  [1].  In*structiire  shock  calculations  were  perfomed  to  nodel  these 
tests,  using  current  sitq>lified  techniques.  Results  were  compared  to  the  test 
data.  Two  (2*D)  finite-elenent  analysis  prograas  were  selected  as  possible 
in- structure  shock  calculation  tools.  These  prograiu  were  also  used  to 
conduct  several  analyses  of  the  test  series.  These  results  were  also  eoapared 
to  the  test  data  and  to  the  siaplified  analysis  results. 

There  were  three  siaplified  analysis  techniques  used  in  this  study.  The 
first  was  Che  seml-eaplrlcal  in-structure  shock  calculation  procedures 
outlined  in  the  US  Aray  IMS  8SS-1  (2].  This  techniqtie  yields  average  in- 
stnicture  shock  motions  for  the  entire  structure.  An  SDOF  analysis  of  the 
wall  closest  to  the  explosion  using  the  US  Aray  Engineer  Vail  Analysis  Code 
(VAC)  [3]  was  the  second  siaplified  technique  examined.  In  this  analysis  the 
aass  and  Che  resistance  of  the  wall  was  replaced  by  an  equivalent  SDOF  system 
and  Che  in- structure  shock  aocions  of  the  wall  were  calculated.  The  third 
siaplified  technique  was  an  SDOF  analysis,  in  which  the  aass  of  Che  entire 
structure  was  lumped  together  at  a  point  and  Che  soil  behind  it  replaced  by  a 
simple  spring.  This  rigid-body  SDOF  analysis  gave  the  overall  aocions  of  Che 
structure  as  a  whole.  These  SDOF  analyses  were  done  on  PC*bised  prograas 
developed  VES., 

The  2-D  finlte-eleaent  analyses  were  accomplished  uslnj;  two  PC-based 
beaa  element  finite -element  analysis  prograas,  STABLE  end  IS!:V3.  STABLE  is  an 
lapllclc  finite -element  program  for  the  dynamic  analysis  of  ixaiaes  subjecteid 
Co  blast  and  ground  shock  loadings  (4] .  This  prograa  is  in  the  public  doaain 
and  was  wriccon  by  JAYCOR,  Vicksburg,  MS,  for  the  US  Amy  Eni;ineer  District, 
Omaha.  1SSV3  is  a  2-D  lumped  parameter  aodel,  explicit,  non], inear,  finite* 
eleiMnc  code  written  by  Applied  Research  Associates,  Southern  Division, 
VicV’sbvirg,  MS  [5].  This  program  was  written  for  VES  in  conjtinction  witl^  this 
project  under  the  direction  of  the  author,  and  is  also  in, Ch«i  public  doaain. 
The  process  of  modeling  this  problea  necessiratad  the  inelusj.on  of  structure- 
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nedla  interaction  (SMI)  for  the  development  of  loads  for  the  striictural 
analysis.  The  development  of  this  SMI  load  was  carefully  evaluated.  Data 
analysis  included  the  examination  of  soil  pressure,  deflection,  velocity,  and 
acceleration  data  and  the  generation  of  shock  spectra. 


1.4  Introductory  Concepts 

The  shock  environment  inside  of  a  dynamically  loaded  structure  is  a 
complex  phenomenon.  The  amplitude,  shape,  and  duration  of  the  acceleration, 
velocity,  and  deflection- time  histories  greatly  Influence  the  response  of 
internal  systems  that  may  be  daisaged  by 
the  shock  environment.  To  simplify  the 
design  and  analysis  of  internal  systems 
against  in- structure  shock,  a  shock 

spectrum  is  generated  to  characterize  MSSS 

the  shock  environment  inside  the  str- 

ucture.  A  shock  spectrum  is  a  plot  of  < 

the  maximum  response  of  an  SDOF  system  ^  S 

in  which  the  base  is  driven  by  the  in- 

structure  shock  environment.  As  can  be  S 

seen  in  Figure  1.1  the  response  of  an  — ' - ^ - - 

SDOF  system  is  dependent  on  the  natural  ^ 

frequency  of  the  system  determined  by  tioil 

the  mass  and  spring  stiffness,  the  damp¬ 
ing  ratio,  and  the  motion  driving  the  Figure  1.1.  Singlo-Degree-of- 

Freedom  System. 

system.  Given  the  motion  and  damping 

parameters  the  maximum  resporice  is  calculated  by  means  of  the  Dxihamel  integral 
[6]  for  a  range  of  natural  frequencies,  and  «  shock  spectra  generated.  For 
many  shock  spbctra  the  natural  frequency  of  the  SDOF  system  is  on  the  abscis¬ 
sa,  while  on  the  ordinate  the  maximum  response,  such  as  peak  spectral  acceler¬ 
ation,  pseudo-velocity ,  or  relative  displacement,  appears.  Peak  spectral 
acceleration  and  pseudo-velocity  are  terms  that  closely  approximate  the  actual 
acceleration  and  velocity  of  the  SDOF  system  while  the  relative  displacement 
is  the  actual  value.  Thus,  fur  a  given  internal  system,  such  as  a  piece  of 
equipment,  all  that  is  needed  to  define  its  maximum  response  to  a  shock 


f  P''se 


Figure  1.1.  Single-Degree-of- 
Freedom  Sys  tern . 
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envlronnent  is  the  natural  frequency  of  the  piece  of  equipment.  All  the 
shock  spectra  generated  in  the  analysis  in  this  study  assumed  no  damping. 

The  shock  spectra,  which  were  generated  in  this  study,  are  a  little  more 
complex  in  that  they  are  plotted  on  tripartite  paper.  There  are  simple  direct 
rdlatioxiships  between  the  maximum  spectral  acceleration,  pseudo-velocity,  and 
relative  displacement  that  allow  these  quantities  to  appear  on  the  same  plot. 
The  log  of  the  natural  frequency  of  the  SDOF  system  appears  on  the  abscissa  of 
these  plots.  There  are  three  ordinates,  one  measuring  the  log  of  psetido- 
velocity  on  the  left  scale,  the  second  the  log  of  the  relative  displacement  on 
a  scale  at  45  degrees*  ^rom  the  vertical  measuring  from  the  lower  right 
comer  to  the  upper  left  comer,  and  the  third,  also  at  45  degrees  from  the 
vertical,  measuring  the  log  of  the  spectral  acceleration  ftom  the  lower  left 
comer  to  the  upper  right  comer.  Numerous  examples  of  shock  spectra  can  be 
found  throughout  this  report.  At  times  in  the  report  spect'.ral  acceleration 
and  pseudo-velocity  are  referred  to  simply  as  acceleration  and  velocity.  The 
prefix  is  implied  whenever  shock  spectra  are  being  discussed. 

In  discussing  shock  spectra,  reference  is  sometimes  made  to  the  constant 
deflection,  velocity  or  acceleration  portion  of  the  plot.  In  the  low  freqtien- 
cies,  the  shock  spectra  of  most  SDOF  systems  tend  to  follow  the  deflection 
scale  progressing  at  a  45  degree  angle  upward  to  th-')  right,  in  the  intermedi¬ 
ate  frequency  region  the  curves  tend  to  be  horizontal  with  a  constant  veloci¬ 
ty,  while  in  higher  frequencies  the  curve  follows  the  acceleration  scale 
downward  to  the  right.  Thc-o  characteristics  are  often  xised  to  generate 
approximate  shock  spectra  with  straight  lines  for  each  portion  of  the  graph. 
Shock  spectra  are  covered  well  in  the  references  [6],  [7],  and  (.8]  and  these 
are  good  references  for  anyone  interested  in  exploring  the  subject  to  greater 
depth. 

Before  going  into  calculation  of  in-structure  shock  and  hence  shock 
spectra,  the  important  concept  of  Hopkinsdn  or  cube- root  scaling  should  be 
introduced.  Host  parameters  in  the  field  of  explosion  effects  are  presented 
as  scaled  by  the  cube  root  of  the  explosive  weight.  This  scaling  relates 
blast  properties  from  an  explosion  of  one  energy  level  to  that  of  an  explosion 

*.\  table  of  factors  for  converting  non- SI  units  of  measurement  to  SI 
(metric)  units  is  presented  on  page  14. 


of  ?  secc.  d  energy  level.  The  energy  is  directly  related  to  the  weight  of  the 
explosive.  This  method  is  widely  covered  in  the  literature  [9].  The  distance 
from  a  bomb  is  scaled  by  dividing  by  the  cube  root  of  the  explosive  weight 
(TNT  equivalent) .  Such  a  scaled  di.«<tance  is  often  referred  to  as  Lambda.  For 
a  given  Lambda,  a  similar  blast  environment  will  result  from  different 
explosive  weights.  Pressures  and  velocities  will  be  identical,  and  time  and 
impulse  are  scaled  by  the  cube  root  of  the  explosive  weight.  Thus,  when  a 
paras ster  is  given  in  a  scaled  form,  it  can  be  applied  to  any  explosive  charge 
size  simply  by  multiplying  by  the  cube  root  of  the  charge  weight.  For  this, 
reason,  most  such  parameters  are  given  in  the  scaled  form.  The  application  of 
cube-root  scaling  was  used  throughout  this  analysis  in  the  development  of 
free- field  loads. 
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CHAPTER  2 
RELATED  TEST  DATA 

2.1  Literature  Search 

An  extensive  literature  search  was  conducted  to  obtain  available  data  to 
assist  in  the  evaluation  of  in- structure  shock  calculations.  It  was  obvious 
that  Buch  of  the  data  available  on  the  subject  of  dynamic  response  of  buried 
concrete  structures  were  not  directly  applicable  to  this  effort.  Most  of  the 
tests  focused  on  plane  wave  type  loadings  and  were  also  largely  concerned  with 
structural  response,  not  in- structure  shock  response.  Also,  most  of  the  data 
was  in  some  way  restricted  from  open  publication.  There  was  one  notable 
exception,  a  series  of  well -controlled  tests  were  conducted  that  included  in¬ 
structure  shock  response  of  b  ^rled  concrete  structures  to  conventional 
explosive  loadings  which  has  been  published  in  the  open  literature.  This 
exception  was  the  CONVEB  test  series. 

2.2  CONVEB  Test  Series 

The  US  Army  Engineer  Va'cerways  Experiment  Station  conducted  the  CONVEB 
tests,  a  series  of  four  backfill  effects  tests  during  the  period  of  March 
through  May  1989.  The  test  procedures  are  covered  extensively  by  Hayes  [1] 
and  [10].  The  analyses  contained  in  these  two  references  were  focused  solely 

I 

on  the  effects  of  backfill  variation  and  structure  media  Interaction  on  the 
structural  response  of  a  buried  test  wall.  In  the  course  of  this  test  series, 
a  large  amount  of  structural  acceleration  data  was  recorded.  These  data, 
along  with  records  of  free-field  acceleration,  free-field  soil  stress,  and 
interface  stress,,  present  a  golden  opportunity  for  the  evaluation  of  in¬ 
structure  shock  calculation  techniques.  For  this  reason,  this  test  series  was 
selected  as  the  primary  focus  of  the  in-structure  shock  analyses  contained  in 
this  report.  The  following  is  a  brief  description  of  the  CONVEB  test  series. 

The  four  tests  in  the  CONVEB  test  series  were  conducted  to  determine  the 
response  of  buried  reinfotced  concrete  walls  to  localized  dynamic  leads  in 
various  backfill  materials  (Table  2.1).  The  first  two  tests  were  conducted  in 
a  clay  backfill  with  two  different  test  walls,  the  second  test  having  a 
stiffer  test  wall  than  the  first.  This  clay  was  characterized  by  a  low  shear 
strength  and  a  low  seismic  Velocity.  The  third  test  was  conducted  with  a 
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Table  2.1.  CONVEB  Backfill  Soil  Propercles. 


1  Test 

Seismic 

Velocity 

(c) 

ft/sec 

Attenuation  I 
Coefficient  1 
(n  or  n) 

CONVEB  1 

122.5 

1100 

2.125 

CONWEB  2 

123.7 

1100 

2.125 

1  CONVEB  3 

116.4 

1000 

3.00C 

1  CONVEB  4 

120 

3000 

1.500 

fl®xible  wa)  1  in  a  sand  backfill  with  a  high  shear  strength  and  low  seismic 
velocity.  The  fourth  test  was  conducted  with  a  flexible  wall  In  a  clay 
backfill  with  low  shear  strength  and  a  high  seismic  velocity.  Figtire  2.1 
through ' Figure  2.4  show  the  test ' structure ,  backfill  properties,  and  the 
properties  of  the  surrounding  ln*situ  soil  for  each  of  the  four  tests. 


Table  2.2. 


CONVEB  Structural 
Material  Properties. 


Table  2 . 1  shows  the  backfill  properties  %rhlch  were  xised  for  all  the  analyses 
in  this  Investigation. 

The  source  of  the  localized  dynam* 
ic  loads  in  each  of  these  tests  was  a 
cylindrical  explosive  charge.  This 
charge  contained  15.4  pounds  of  C4  ex¬ 
plosives  in  a  closed  steel  case  27  inch¬ 
es  long  with  an  inside  diameter  of  3.548 
inches  and  a  thickness  of  0.166  inch. 

The  charge  was  oriented  vertically  in 
CONVEB  1,  3,  and  4,  and  horizontally  In 
CONVEB  2. 

The  test  article  consisted  of  a 
reusable  reaction  structure  supporting 
the  test  wall.  The  test  wall  was  bolted 
to  the  heavily  reinforced  concrete  reac¬ 
tion  structure,  forming  a  relatively 

rigid  joint.  Dimensions  and  gage  layout  of  the  reaction  structure  are  shown 


Test 

Component 

f’ 

(psi) 

Steel 

Yield 

(psi) 

CONVEB  1 
Test  Vali 

6095 

67424 

CONVEB  2 
Test  Vail 

6398 

67424 

CONVEB  3 
CONVEB  4 
Test  Valla 

5855 

67424 

Reaction 

Structure 

6398 

67424 
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in  F.lgure  2.5.  Test  walls  for  CONWEB  1,  3.  and  4  had  a  length- to -thickness 
<L/t)  ratio  of  10  while  the  L/t  ratio  of  the  wall  in  CONWEB  2  was  5.  The 
designs  of  the  test  walls  are  shown  in  Figure  2.6  and  Figure  2.7,  and  the  test 
wall  gage  layout  is  shown  in  Figure  2.8.  Material  properties  are  presented  in, 
Table  2.2, 


Figure  2.1.  Section  Through  CONWEB  1  Test  Bed,  Showing  Backfill 
Properties. 
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Figure  2.3.  Section  Through  CONVEB  3  Test  Bed,  Shoving  Backfill 
Properties . 


Figure  2.4.  Section  Through  CONUEB  4  Test  Bed,  Shoving  ln*Situ  Soil 
Properties . 
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65.19 


Figure  2.5.  Reaction  Structure  Dimensions  and  Gage  Layout. 


25 


Figure  2.7.  CONWEB 
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Figure  2.8‘.  CONVEB  Test  Vail  Instrunentation. 
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CHAPTER  3 

SIMPLIFIED  ANALYSIS  TEaiNIQUES 


3.1  TM5  855-1  In-Structixre  Shock  Procedure 

A  current  method  for  the  calculation  of  in- structure  shock  is  the  semi- 
empirical  procedure  in  the  Amy  Technical  Manual  (TM)  5  855-1  [2].  The 
procedure  begins  with  the  calculation  of  the  average  free-field  accelerations, 
velocities,  and  deflections.  Correction  factors  are  then  applied  to  convert 
these  free-field  values  to  in- structure  shock  values.  The  following  is  a 
brief  explanation  of  the  method.  Later,  the  procedure  will  be  applied  tc  the 
CONVEB  test  ser5.es. 

The  most  important  parameters  controlling  the  free-field  motion  are:  the 
source  of  the  ground  shock,  the  mechanical  properties  of  the  free-field,  and 
the  range  (distance)  to  the  point  of  interest.  The  sottrce  of  the  grotind  shock 
is  detemlned  by  the  size  and  depth  of  burial  of  the  botd>. 

There  are  many  types  of  explosives  used  in  bombs,  each  having  vary 
different  characteristics.  In  order  to  simplify  analysis,  the  weight  of 
explosive  is  converted  to  an  equivalent  TNT  weight  using  equivalency  factors. 
Unfortunately,  these  conversion  factors  were  originally  generated  for  free  air 
(above  ground)  explosions  and  in  low  pressure  ranges Application  of  these 
factors  to  belowground  explosions  is  therefore  open  to  criticism.  This  is, 
however,  common  practice  and  is  the  best  available  approach. 

After  conversion  to  an  equivalent  TNT  charge  weight,  the  affect  of  depth 
of  burial  must  also  be  taken  into  account.  It  is  obvious  that  a  bomb,  which 
is  barely  penetrating  into  the  soil,  will  not  be  as  effective  in  transmitting 
energy  into  ground  shock  as  would  a  bomb  that  is  fully  buried.  A  ground  shock 
coupling  factor  is  used  to  account  for  this  phenomenon.  This  single  factor  is 
applied  to  all  ground  shock  parameters  and  is  simply  the  ratio  of  the  shallow¬ 
er  buried  bomb's  parameters  to  that  of  a  fully  burled  bomb,  as  shown  in 
Equation  (1).  The  shock  parameters  in  question  are:  pressure  (P),  velocity 
(V),  deflection  (d) ,  impulse  (I),  and  acceleration  (a). 


f- 


{P,  V,  d,  I,  a)  near  surface 
(P.  V,  d,  I,  a)  contained 


(1) 
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The  propagation  of  shock  through  soil  is  a  conplex  problaa,  dependent  on 
a  largo  nunber  of  variables.  At  the  present  time,  the  TM5  855-1  procedure  for 
the  calculation  of  free -field  nwtlons  contains  only  two  explicit  variables  for 
the  characterization  of  the  soil.  These  variables  are  the  seismic  velocity 
(c)  and  the  attenuation  coefficient  (n).  Seismic  velocity,  c,  is  used  as  ^n 
Index  of  soil  properties  for  ground  shock  prediction,  providing  an  indication 
of  the  soil  stiffness  (M)  and  mass  density  (p),  as  shown  in  Equation  (2).  The 

(2) 

value  of  c  is  also  strongly  dependant  on  the  degree  of  saturation  of  a 
cohesive  soil.  In  general,  the  higher  the  val\ie  of  c,  the  better  the  soil  is 
at  transmitting  shock.  The  reverse  is  true  for  the  attenuation  coefficient, 
n,  .  %diich  is  a  measure  of  the  energy  used  in  the  irreversible  cmshing  of  the 
soil  voids.  A  high  value  of  n  indicates  a  soil  that  will  quickly  attenuate  a 
shock  with  distance.  The  characterization  of  a  soil  with  only  two  parameters 
is  an  extreme  simplification  and  should  be  used  with  great  caution. 

a»f’.50fc,(^)‘*‘*‘*  (3) 


1 


(5) 


The  final  parameter  used  in  the  TH5  855-1  procedure  for  the  calculation 
of  free- field  motion  is  the  range  (R)  or  distance  of  the  bomb  from  the  point 
of  interest.  As  would  be  expected,  ground  shock  and  free -field  motions 
decrease  with  increasing  range.  The  free-field  peak  acceleration  (a)  in  g's, 
peak  particle  velocity  (V)  in  ft/sec,  and  peak  displacement  <d)  in  feet  can 
now  be  calculated  using  Equations  (3),  (4)  and  (5).  The  weapon  yield  {U)  is 
in  pounds  TNT  equivalent,  coupling  factor  (f)  and  attenuation  coefficient  (n) 
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as  defined  above,  seismic  velocity  (c)  in  ft/sec,  and  range  (R)  in  feet.  Note 
that  Hopkinson's,  cube -root  scaling  was  used  extensively  in  the  development  of 
these  eqviations . 


Figure  3.1.  Average  Free-Field  Motion  Configuration  for  Side  Burst 
Load  Case . 

This  brings  us  now  to  the  actual  calculation  of  in-structure  shock. 
Figure  3.1  shows  the  general  configuration  of  the  side  burst  loading  of  a 
rectangular  structure.  The  two  ranges,  and  R2,  are  the  distance  in  feet 
from  the  front  and  the  back  of  the  structure,  respectively.  These  two  ranges, 
along  with  the  parameters  discussed  above,  are  all  that  is  required  for  the 
calculation  of  the  average  acceleration  across  the  structure.  Equations  (6), 
(7),  and  (8)  are  used  to  calculate  the  average  acceleration,  velocity,  and 
deflection  across  the  structure.  These  equations  were  developed  by  integrat¬ 
ing  the  acceleration- range  relationship  in  Equation  (3)  across ’  the  stnicture 
and  finding  a  uniform  acceleration  that  yields  the  same  integral  across  this 
range .  , 

.  (6) 


31 


(7) 


(a-l)  (Rt-Rr) 

(n-H 

500ir¥  (V**»-J!i-***)  (Z) 

“  c(n-2)  (/2,-Jii) 

Given  the  average  free -field  awtions  calculated  above,  the  next  step  is 
to  modify  these  values,  applying  empirical  factors  to  produce  an  approximation 
of  in- structure  shock  motions.  A  reduction  factor  (RF)  is  calctxlated  based  on 
the  geometry  of  the  structure,  vail  height  or  width,  structure  length,  and 
range.  A  nomograph  is  \ised  to  find  this  empirical  relationship.  This 
reduction  factor  is  applied  to  accelerations  and  velocities.  These  values  are 
taken  as  an  approximation  of  the  accelerations,  velocities,  and  deflections, 
which  would  be  measured  inside  the  structure  during  an  event  under  these 
conditions. 

The  final  step  in  the  TM5-855-1  in-structure  shock  analysis  procediire  is 
the  generation  of  a  shock  spectrum.  As  discussed  above,  a  shock  spectrum  is  a 
very  useful  tool  for  the  analysis  of  the  dynamic  response  of  equipment  and 
other  subsystems  inside  a  structure.  Generation  of  shock  spectrum  is  simply  a 
matter  of  multiplying  the  above  in- structure  motions  by  empirical  amplifica¬ 
tion  factors.  These  factors  are  1.2,  1.5,  and  1.6  for  displacement,  velocity, 
and  accelerations,  respectively.  Response  values  are  plotted  on  the  tripar¬ 
tite  shock  spectrum  curves  as  the  maximum  relative  displacement,  maxlimni 
pseudovelocity,  and  maximum  acceleration.  These  shock  spectra  can  be  directly 
compared  to  information  on  the  response  of  equipment  inside  the  structure. 
Equipment  response  is  often  presented  es  a  shock  spectrum  which  gives  recom¬ 
mended  limits  on  equipment  motion.  Overlaying  the  structural  shock  spectrtim 
on  these  equipment  fragility  curves  allows  an  estlMte  of  expected  damage  and 
the  need  for  shock  isolation. 

3.2  TMS  855-1  CONVEB  Calculations 

Table  3.1  shows  the  values  used  for  the  calculation  of  average  free- 
fleld  notions  for  CONVEB  tests  1,  2,  3,  and  4.  These  values  were  input  into 
Equations  (6),  (7),  and  (8).  Table  3.2  shows  the  average  in- Structure  notions 
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Table  3.1.  CONWEB  TM5  855-1  Free-Field  Response  Parameters. 


T'  st 

V 

lbs. 

f 

n 

c 

ft/sec 

Ri 

ft 

R2 

ft 

RF 

CONVEB  1 

19.71 

Bsm 

1100 

5.0 

9.717 

.48 

CONVEB  2 

19.71 

mm 

2.25 

1000 

5.0 

9.717 

.48 

CONVEB  3 

19.17, 

1.0 

3.0 

1000 

5.0 

9.717 

.48 

CONVEB  4 

19.71 

1.0 

1.5 

3000 

5.0 

9.717 

.48 

and  spectral  values  for  each  test.  The  shock  spectrum  values  are  shown 
plotted  against  shock  spectra  generated  from  test  data  in  Figure  3.2, 

Figure  3.3,  Figure  3.4,  and  Figure  3.5. 

The  TMS  855-1  procedure  yielded  very  good  results  for  CONVEB  1  and  2  and 

Table  3.2.  CONVEB  Tl^  855-1  Average  In-Structure  Response  and  Spectral 
Values . 


Test 

A 

g 

V 

in. /sec 

D 

in. 

1.6A 

g 

1.5V 
in. /sec 

1.2D 

in. 

CONVEB  1 

539 

114 

4.98 

862 

170.0 

5.98 

CONVEB  2 

439 

110 

1.80 

702 

166 

2.16 

1  CONVEB  3 

614 

110 

1.29 

983 

166 

1.54  1 

1  CONVEB  4 

2549 

219 

3.32 

4078 

328 

3.98  1 

gave  a  conservative  estimate  of  response  for  CONVEB  3.  The  relatively  good 
results  of  these  calculations  are  not  surprising,  as  the  procedure  was 
developed  from  testa  on  simple  box  structures  similar  to  the  CONVEB  test 
articles.  Figure  3.5  shows  that  CONVEB  4  was  undetpredicted  according  to  gage 
ASHIO  and  overpredicted  for  gage  ASHll.  If  gage  ASHIO  is  correct  and  CONVEB  4 
was  underpredicted,  this  could  demonstrate  a  problem  often  seen  in  dynamic 
structural  analysis,  the  problem  of  calculating  the  response  of  structures  in 
high  seismic  velocity  clay.  However,  another  explanation  could  be  the 


variability  in  the  testing  itself.  CONWEB  4  was  an  in-sicu  test  in  which  the 
soil  was  excavated  and  the  structure  was  placed  in  the  excavation.  Difficulty 
was  encountered  in  assuring  a  proper  placement  of  backfill  around  the  struc¬ 
ture  to  fill  the  gap  between  the  structure  and  the  in-situ  soil.  There  was 
also  uncertainty  in  the  level  of  the  ground  water  table  in  this  test.  Both  of 
these  variables  could  drastically  affect  the  response  of  the  structure. 
Improper  backfill  would  allow  unrestrained  Mtion  of  the  strxicture,  and  the 
presence  of  water  could  greatly  Increase  the  dynamic  load.  Gage  failure  may 
be  a  better  explanation  of  the  apparent  underprediction.  If  gage  ASHll  is 
correct,  the  IMS  RS5-1  gave  a  conservative  estimate  of  horizontal  in- structure 
shock. 

The  semi -empirical  calculational  procedure  in  TM5  855-1  has  now  been 
demonstrated  to  be  capable  of  giving  a  reasonable  estimate  of  the  in-structtire 
shock  response  of  the  iiist  three  CONVEB  tests.  Inconsistences  in  the  CONVEB 
4  test  results  made  it  difficult  to  evaluate  the  analysis  results  for  this 
test.  The  test  was  underpredicted  when  compared  to  one  gage,  while  the 
aralysls  gave  reasonable  results  when  compared  to  another  gage  at  almost  the 
same  location.  This  test  nay  have  been  overloaded  and/or  underrestrained  due 
to  testing  problems  or  to  have  had  gage  problems. 

The  greatest  weakness  of  this  procedure  is  that  it  is  sn  empirical 
method  developed  for  simple  box-like  structures.  The  assumptions  that  must  be 
made  to  apply  this  procedure  to  more  complicated  structures  have  not  l!>een 
thoroughly  evaluated.  Great  caution  should  be  used  in  the  application  of  such 
an  empiritally  derived  method  to  sitxiations  outside  of  those  sctually  tested. 
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Figure  3.2.  CONWEB  1,  TM5  855-1  Generated  Shock  Spectmai  ve.  Teet  Data. 
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Figure  3.3.  CONVEB  2,  TM5  855-1  Generated  Shock  Spectrua  vs.  Test  Data. 
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Figvire  3.4.  CONUEB  3,  TM5  8S5>1  Geiuiratfld  Shock  Spectrua  vs.  Test 
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Figure  3.5.  CONVEB  4,  TMS  855-1  Cenersted  Shock  Spectrua  vs.  Test 


3 . 3  SDOii'  Analyses 

Tne  single -degree -of- freedom  (SDOF)  analysis  Is  the  second  simplified 
method  examined  for  in-str’acture  shock  calculations.  This  method  consists  of 
the  reduction  of  the  problem  to  a  simple  spring,  mass,  and  damper  system. 
Application  of  this  type  of  analysis  is  widely ' covered  in  the  literature  [1], 
[7].  In  an  SDOF  model,  the  spring  Is  a  resistance  element  which  models  the 
static  resistance  of  the  structural  element  of  interest.  The  mass  and  damper 
ase  selected  so  that  the  resulting  system  will  have  the  same  frequency  and 
damping  characteristics  as  the  prototype  structxire. 

The  SDOF  analysis  procedure  was  used  to  model  the  CONVEB  tests  in  two 
ways,  the  wall  facing  the  bomb  was  analyzed  as  one  SDOF  system,  and  the 
horizontal  rigid-body  motion  of  the  entire  structure  was  modeled  in  a  second 
decoupled  SDOF  analysis.  Output  frem  the  wall  analysis  was  compared  to  data 
from  gages  on  the  wall ,  while  the  rigid-body  analysis  results  were  compared  to 
data  from  internal  gages  on  the  floor.  This  Ignores  the  contribution  of  the 
rigid-body  motion  to  the  response  of  the  wall.  Since  the  peak  wall  response 
occurs  well  before  the  peak  rigid-body  response,  neglecting  of  the  rigid-body 
contribution  is  thought  to  be  a  reasonable  assumption. 

3.3.1  Front  Wall  SDOF  Analysis 

The  front  wall  SDOF  aiialysls  was  carried  out  using  the  Corps  of  Engi¬ 
neers  PC-based  Wall  Analysis  Code  (WAC)  [3].  WAC  uses  the  procedures  in  TM5 
855-1  to  develop  the  load  mass  factors  for  a  given  wall.  Hultilinear  resis¬ 
tance  functions  are  computed  based  on  yield- line  theory.  The  resistance 
functions  for  the  test  walls  in  the  CONWEB  tests  are  shown  in  Figure  3.6. 

Note  that  the  calculated  resistance  functions  were  the  same  for  CONWEB  1,  3, 
and  4,  as  the  wall  designs  were  the  sane  for  these  tests,  and  the  same 
material  properties  were  used. 

The  wall  loading  used,  was  a  modification  of  the  free-field  stresses 
calculated  by  the  procedures  in  TM5,  855-1  (2].  The  soil  properties  used  are 
shown  in  Table  2.1  (in  chapter  11  of  this  report).  Free-field  stresses  were 
modified  by  applying  a  reflected  pressure  factor  of  1.5  to  the  beginning  of 
the  pressure- time  history.  This  models  the  buildup  of  reflected  pressure  on 
the  structure  before  a  tensile  relief  wave  can  arrive  from  the  edge  of  the 
structure,  or  before  wall  movement  can  relieve  the  stress.  According  to  TM5 


37 


855-1  the  duration  of  this  reflected  pressure  for  these  tests  should  be  taken 
as  six  wave  transit  times  through  the  thickness  of  the  wall,  0.2  msec  in 
CONWEB  1,  and  0.4  msec  in  CONWEB  2,  3,  and  4.  Applying  this  reflected 
pressure- time  history  over  the  entire  surface  of  the  test  wall  is  overly 
conservative,  due  to  the  nonuniform  nature  of  conventional  explosive  loadings. 
The  entire  pressure- time  history  was  reduced  by  applying  a  uniform  load  factor 
of  0.71  as  per  TH5  855-1,  calculated  based  on  the  aspect  ratio  of  the  wall  and 
the  weapon  range.  The  final  pressure- time  histories  used  in  the  analysis  are 
shown  in  Figure  3.7  through  Figure  3.10. 

Peak  deflections,  velocities,  and  accelerations  resulting  from  the  wall 

Table  3.3.  Peak  Response  Output  from  SDOF  Vail  Analysis. 


Test 

Deflection 

(in.) 

Velocity 
(in. /sec) 

Acceleration  | 

(g) 

CONVEB  1 

60 

3050 

2095  1 

CONVEB  2 

7.1 

828 

874.0  1 

CONVEB  3 

0.68 

187 

443.2  1 

CONVEB  4 

25 

2310 

3391  1 

SDOF  analysis  are  shown  in  Table  3.3.  As  in  the  TM5  855-1  analysis  above,  a 
factor  of  1.2,  1.5,  and  1.6  was  applied  to  the  peak  deflection,  velocity,  and 
acceleration  respectively  to  create  shock  spectrum  values.  The  shock 

spectra  generated  are  shown  in  Figure  3.11  through  Figure  3.14.  The  SDOF 
procedure  overpredicted  the  peak  deflection  and  velocity  portions  of  the  shock 
spectra  and  correctly  predicted  the  acceleration  portion  for  CONVEB  1.  CONVEB 
2  results  show  a  very  reasonable  comparison  to  test  data,  while  the  CONVEB  3 
analysis  shows  an  underprediction  of  the  spectral  values.  CONVEB  4  was 
reasonably  well  predicted  for  all  three  parts  of  the  spectmm. 

The  results  of  the  wall  SDOF  analysis  are  relatively  inconsistent  and 
reflect  a  problem  often  encountered  in  this  type  of  work;  that  is,  the 
difficulty  of  developing  reasonable  loads  to  apply  to  the  SDOF  model.  In  work 

by  Hayes  [1],  In  which  a  similar  SDOF  structural  analysis  of  the  CONVEB  test 

\ 
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vails  gave  similar  anomalous  results,  a  major  conclusion  was  that  seismic 
velocity  alone  is  not  sufficient  to  characterize  a  given  backfill  material. 
The  loads  generated  by  the  TM5  855-1  procedures,  which  were  used  here,  are 
very  dependent  bn  seismic  velocity.  Also,  it  is  thought  that  the  TH5  855-1 
procedures  do  not  adequately  model  structure -media  interaction;  that  is,  the 
reflection  of  the  load  at  the  soil -structure  interface  and  the  relative 
movement  of  the  soil  and  the  structure.  There  is  current  research  at  WES  to 
address  these  deficiencies. 
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Figure  3.12.  CONWEB  2,  SDOF  Wall  Analysis  Shock  Spectrua  vs.  Test  Data. 
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Figure  3.13.  CONWEB  3,  SDOF  Wall  Analysis  Shock  Spectnia  vs.  Test  Data. 
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Flgiire  3.14.  CONVEB  4.  SDOF  V«ll  Analysis  Shock  Spectrua  vs.  Test  Data. 


3.3.2  Rigid-Body  SDOF  Analysis 

A  simple  approximation  of  the  horizontal  shock  environment  on  the  floor 
of  a  structure  is  the  overall  rigid-body  motion  of  the  entire  structure. 
Rigid-body  motion  was  analyzed  by  reducing  the  problem  to  an  SDOF  system,  with 
the  total  mass  of  the  structure  concentrated  to  a  single  point,  and  with  the 
soil  behind  the  structure  acting  as  a  simple  linear  spring.  The  load  used  was 
the  same  as  that  calculated  for  the  wall  SDOF  above.  Shock  spectra  were 
generated  from  the  resulting  motion  histories  and  compared  to  those  from 
internal  acceleration  gages  located  on  the  structxire's  floor.  It  should  be 
noted  that  this  SDOF  model  only  gives  meaningful  results  up  to  the  initial 
peak  positive  deflection.  The  procedure  does  not  model  the  complex  interac¬ 
tions  which  take  place  at  later  times  as  the  structure  rebounds.  However,  the 
highest  deflections,  velocities,  and  accelerations  almost  always  take  place 
before  this  happens  and  this  limitation  does  not  hamper  the  generation  of 
shock  spectra. 

The  soil -spring  constant  was  calculated  using  a  procedure  presented  by 
Whitman  and  Richart  [11]  for  dynamically  loaded  foundations.  Equation  (9)  was 
used  to  compute  the  soil-spring  constant  for  elastic  response  of  the  soil  at 
the  back  of  the  structure,  where  G  is  the  soil's  shear  modulus  (psi);  B,  a 

(9) 

dimensionless  aspect  ratio  coefficient;  p  the  Poisson's  ratio  for  the  soil; 
and  B  and  L  the  height  and  width  of  the  structure  respectively  (inches).  The 
shear  modulus  G,  was  calculated  from  the  soil  density  (p)  and  shear  wave 
velocity,  V,,  using  Equation  (10).  The  shear  wave  velocity  was  calculated 
using  Equation  (11)  from  the  compressive  seismic  velocity.  Table  3.4  shows 
the  valites  used  for  this  calculation  and  the  resulting  soil-spring  constant; 
the  total  mass  of  the  structure  is  also  shotm.  The  Poisson’s  Ratios  used  were 

G  -  pV,*  (10) 

yl  .  y3  (11) 

•  '  2(l-p) 

chosen  based  on  soli  tests  conducted  in  conjunction  with  the  CONVEB  test 
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saries  [1].  These  values  are  relatively  hi^  but  are  within  the  range  of 
those  suggested  by  Bowels  [12] •  Frictional  forces  that  exist  at  the  top  and 
bottom  of  the  structure  are  neglected  in  this  analysis.  A  relatively  low 
damping  value  of  one  percent  of  critical  was  included  for  stability  and  does 
not  affect  the  first  peak  response. 

Given  the  above  mass,  soil* spring  constant,  damping  coefficient,  and  the 
load  calculated  earlier,  the  rigid-body  SDOF  analysis  was  completed  iising  a 
computerized  version  of  the  procedures  developed  by  Biggs  [7].  Velocity-time 
histories  up  to  the  time  of  positive  deflection  were  input  into  the 

same  shuck- spectra  generation  program  used  for  the  test  data.  Comparisons  of 
these  SDOF  spectra  with  spectra  generated  from  test  data  are  shotm  in 
Figure  3.15  through  Figure  3.18.  As  can  be  seen  in  these  figures,  the  rigid- 
body  SDOF  analysis  gave  reasonable  predictions  of  peak  deflections  and 
accelerations  for  CONVEB  1,  2,  and  3.  Peak  velocities  were  overpredicted  for 
these  saoKJ  tests.  The  same  trends  can  be  detected  in  the  CONVEB  4  predictions 
in  comparison  with  gage  ASHll.  In  comparison  to  ASHIO  the  CONVEB  4  predic¬ 
tions  were  low  for  deflection  and  acceleration  and  reasonable  for  velocity. 

Table  3.4.  Rigid- Body  SOOF,  Soil -Spring  Parameters. 


Test 

■ 

n 

G 

psi 

B. 

k 

Ib/in. 

M 

lb -sec*/ In. 

CONVEB  1 

.45 

332 

2906 

2.22 

1.269xl0‘ 

73.21 

CONVEB  2 

.45 

332 

2935, 

2.22 

1.281x10* 

84.48 

CONVEB  3 

.33 

504 

6370 

2.22 

2.283x10* 

73.21 

CONVEB  4 

.45 

904 

21176 

2.22 

9.245x10* 

73.21 

indicating  that  the  overall  response  in  this  test  was  underpredicCed.  Under¬ 
prediction  of  CONVEB  4  in  comparison  with  gage  ASHIO  has  occurred  in  each  of 
the  analysis  procedures  used  so  far,  perhaps  lending  more  credence  to  the 
possibility  of  testing  or  gage  problems  as  discussed  above. 

The  consistent  overprediction  of  peak  spectral  velocity  in  the  rigid- 
body  SDOF  analysis  is  more  difficult  to  explain.  Possible  overloading  of  the 
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structure  due  to  problem  Ip  modeling  structure  media  interaction,  as  dis¬ 
cussed  above,  has  more  effect  on  the  total  impulse  imparted  into  structure 
than  it  does  on  the  peak  pressure  experienced  at  the  front  vail.  The  peak 
acceleration  of  an  SDOF  model  is  very  dependent  on  peak  pressure,  while  the 
peak  velocity  is  more  dependent  on  the  impulse.  Thus,  overloading  the 
structure  with  Impulse  leads  to  overprediction  of  velocity.  This  should  also 
lead  to  an  overptediction  of  deflection.  The  fsct  that  deflection  predictions 
are  reasonable,  indicates  a  compensating  overrestraint  of  die  model,  perhaps 
because  the  Poisson's  Ratios  assumed  in  the  analysis  were  on  the  high  pnd  of 
expected  values.  Taken  together,  all  this  illustrates  the  difficulty  of 
analyzing  the  complex  response  of  a  buried  structure  using  a  simple  SDOF 
model.  It  should  be  noted,  however,  that  the  velocity  predictions  were 
conservative. 
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ire  3. IS.  CONVEB  1,  Rlgld*Body  SDOF  Analysis  Shock  Spoctrua  vs.  Test 
Data. 
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Figure  3.16.  CONVEB  2,  Rigid- Body  SDOF  Analysis  Shock  Speetrua  vs.  Test 
Data. 
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Figure  3.1/.  COHVEB  3,  Rigid-Body  SDOF  Analysis  Shock  Speecrua  vs.  Test 
Data. 
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Figure  3.18. 
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CHAPTER  4 

IMPLICIT  FINITE- ELEMENT  ANALYSIS  TECHNIQUE 

4.1  STABLE  Program 

STABLE  is  an  implicit  finite-element  program  for  the  dynamic  analysis  of 
frames  subjected  to  blast  and  ground  shock  loadings.  This  {Program  is  in  the 
public  domain  and  was  written  by  JAYCOR.  Vicksburg,  MS,  for  the  US  Amy 
Engineer  District,  Omaha.  STABLE  was  in  existence  at  the  beginning  of  this 
project  ard,  as  a  validated  program,  was  selected  as  a  candidate  fcr  incorpo¬ 
ration  into  an  in-structure  shock  design  tool.  This  program  was  not  specifi¬ 
cally  designed  for  in- structure  shock  calculations,  and  the  following  is  an 
evaluation  of  the  application  of  STABLE  to  this  class  of  problems. 

4.2  STABLE  Program  Formulation 

Extensive  documentation  of  the  formulation  and  validation  of  the 
computer  code,  STABLE,  has  been  provided  by  Bryant,  Campbell,  Smith,  and 
Flathau  (4],  (13),  and  (14].  The  following  is  a  brief  overview  of  the 
formulation  of  STABLE  as  it  applies  to  the  evaluation  of  the  program  as  an  in- 
structure  shock  tool. 

The  general  formulation  of  STABLE  is  an  Implicit  finite -element  analysis 
program.  Implicit  formulation  refers  to  the  solution  method  used  to  solve  the 
equations  of  motion  as  the  problem  proceeds  through  time.  In  this  case  the 
Ne%niark  integration  scheme  which  is  widely  covered  in  the  literature  [16]  was 
used.  As  such,  the  implicit  solution  of  the  structural  equations  is  xmcondi- 
tionally  stable;  that  is,  the  choice  of  a  large  time  step  will  ^t  cause  the 
solution  to  go  unstable.  This  large  time  step  is  the  main  advantage  of  such  a 
formulation.  However,  with  each  time  step  the  solution  of  the  equations  of 
motion  requires  a  relatively  large  amount  of  calculational  effort  and,  hence, 
computer  time.  An  alternative  formulation  of  a  finite -element  analysis 
program  is  an  explicit  type,  in  which  the  time  step  must  be  small  to  assure 
stability  but  the  computational  effort  is  smaller  at  each  step.  The  program, 
ISSV3,  which  will  be  examined  later,  is  of  an  explicit  formulation. 

The  elements  available  in  STABLE  are  one -dimensional,  prismatic,  beam 
elements  of  constant  cross  section,  moment  of  inertia,  and  plastic  moment 
capacity.  The  material  model  available  in  STABLE  is  a  general  bilinear, 
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elasCO'pIastic  material.  A  steel  shape  table  allows  direct  input  of  standard 
steel  shape  properties.  Reinforced  concrete  is  modeled  by  use  of  an  equiva¬ 
lent  steel  section,  inputting  equivalent  areas,  moment  of  inertia,  and  plastic 
moment  capacity.  A  reinforced  concrete  material  model,  which  appears  in  the 
STABLE  manual,  has  implementation  problems  that  did  not  allow  its  use. 

Dynamic  loading  models  used  in  this  analysis  Included  direct  input 
pressure- time  histories  and  loads  developed  through  structure  media  interac¬ 
tion  (SMI).  Pressure* time  histories  that  were  input  were  simplified  waveforms 
taken  from  the  interface  pressure  gage  data  at  the  front  wall  of  the  stzoic- 
ture.  Restraint  of  the  horizontal  motion  of  the  structure  was  provided  by; 
frictional  forces  modeled  by  Coulomb  dampers  at  the  top  and  bottom  of  the 
structure,  and  SMI  loads  at  the  back  or  the  structure.  SMI  was  also  modeled 
at  the  top  and  bottom  of  the  structure  for  vertical  reactions.  STABLE 
calculates  SMI  loads  using  a  simplified  approach  as  shown  in  Equations  (12) 
and  (13). 

±  pCi(V„„-V„)  o„>0 

a„  ■  SMI  normal  stress. 

p  ■  Soil  mass  density. 

•  Praia- field  compression  wave  speed, 

■  Normal  free- field  soil  velocity. 

V„  •  Normal  structure  velocity.  . 

|t|  i  c+o„tan^ 

T  ■  SMI  shear  stress. 
p  ■  Soil  mass  density, 
c,  •  Free- field  shear  wave  speed, 

Vfft  ■  Teuigential  free- field  soil  velocity. 

V,,  ■  Tangential  structure  velocity. 

^  ■  Soil  friction  angle. 

Normal  and  shear  stresses  resulting  from  the  SMI  model  are  a  combination  of 
the  free- field  stresses  and  the  stresses  due  to  the  relative  movement  of  the 
free- field  and  the  structure.  In  this  analysis,  no  free -field  stresses  were 
applied  at  the  nodes  where  SMI  was  being  calculated.  Thus,  only  the  relative 
motion  stresses  were  in  existence  at  these  points.  The  application  of  actual 
interface  loads  at  the  front  face  already  includes  all  SMI  effects  in  the 
data. 


(12) 


(13) 
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4.3  STABLE  CONVEB  Analysis 

The  evaluation  of  STABLE  as  an  in-structure  shock  analysis  technique 
focused  on  the  analysis  of  the  CONUEB  2  test.  It  was  planned  that  if  this 
evaluation  was  favorable,  the  other  CONVEB  tests  would  be  analyzed.  As  will 
be  discussed  below,  the  results  of  CONVEB  2  analysis  led  to  the  rejection  of 
this  program  as  an  in- structure  shock  analysis  tool,  and  no  further  analyses 
were  conducted.  The  analyses,  which  were  conducted  on  the  CONVEB  2  test, 
included  a  coarse  grid  analysis,  a  fine  grid  analysis,  and  a  fine  grid 
analysis  with  soil  springs  replacing  the  SMI  at  the  back  of  the  structure. 

The  following  is  a  discussion  of  the  results  of  each  of  these  analyses: 

4.3.1  CONVEB  2  Coarse  Grid  Analysis 

For  this  analysis,  a  relatively  coarse  two-dimensional  finite-element 
grid  was  generated  to  model  CONVEB  2  (Figure  4.1).  The  model  consists  of  a  1 
foot- thick  slice  taken  through  the  centerline  of  the  structure.  All  major 
reinforcement  in  the  structure  are  in  this  plane,  and  the  effects  of  all  out- 
of-plane  reinforcement  was  neglected. 

Equivalent  steel  sections  were  used  to  model  the  concrete  cross  sec¬ 
tions.  Equation  (14)  was  used  to  find  the  equivalent  area  and  is  taken  from 
TM5-1300  [16). 

A  m  Equivalent  Area  (inches^) 

Eg  ■  Steel  Youngs  Modulus  {psi)  (14) 

Eg  ■  Concrete  Youngs  Modulus  {psi) 

A,  ■  Steel  Area  {inches*) 

Ag  ■  Gross  Concrete  Area  {inches*} 

Equation  (IS)  was  used  to  find  the  thrust  capacity  and  came  from  the  sane 
source  [16). 

Pg  •  0.9Sfd\{Ag-A,)->’A,fd, 

Pg  •  Compressive  Capacity  {psi) 

fd' g  •  Dynamic  Concrete  Strength  {psi}  (15) 

fdg  •  Dynamic  Steel  Yield  Strength  {psi) 

Ag  m  Gross  Area  {inches*) 

A,  •  Steel  Area  {inches*) 

Equivalent  moment  capacity  is  given  by  Equation  (16)  also  from  IMS  1300  [16]. 
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(16) 


if,  -  (A,~A\)fd,(d-a/2)*A*^d,{d-d') 
a  -  (A,-A,)/d^(0.85Jbfd„) 

>  Moment  Capacity  {in. -lb) 
fd,  ■  Dynamic  Steel  Yield  Strength  ipai) 
fdg  •  Dynamic  Concrete  Strength  ipsi) 
b  ■  Section  Width  {inches) 

A',  "  Tensile  Steel  Area  (inches*) 

A,  ■  Compressive  Steel  Area  (inches*) 


Equation  (17),  Biggs  [7],  gives  the  equivalent  aonent  of  inertia. 

I  -  -^(s.sp-^o.osa) 

I  ■  Moment  of  Inertia  (inches*)  , 
be  Section  Width  (inches) 
d  ■  Depth  to  Tension  Steel  (inches) 
p  ■  Reinforcement  Ratio 


The  resulting  equivalent  section  properties  are  shovrh  in  Table  4.1. 

Table  4.1.  Section  Properties,  CONUEB  2  STABLE  Analysis. 


1  Section 

Compressive 

Capacity 

lb 

Moment 
Capacity 
in. -lb 

Moment  of 

Inertia 

inches* 

Equivalent  I 
Area  1 
inches*  I 

•  Front 

H  Wall 

677358 

224515 

303 

109.2 

1  Reaction 

1  Structure 

844401 

603711 

870 

144 

,  The  results  of  the  coarse  grid  analysis  are  presented  in  Figure  4.2 
through  Figure  4.9.  Figure  4.2  shows  a  comparison  of  the  horizontal  accelera¬ 
tion  record  produced  by  STABLE  at  node  10,  compared  to  the  test  data  at  that 
same  location.  The  initial  peak  values  of  acceleration  coispare  reasonably 
well.  However,  a  late-time,  high-frequency  oscillation  occurs  in  the  STABLE 
output,  which  is  not  seen  in  the  test  data.  This  oscillation  or  noise  in  the 
acceleration  output  is  disturbing,  but  not  fatal,  in  an  in-structure  shock 
analysis.  The  integration  of  the  acceleration  record  to  yield  the  velocity 
record  (Figure  4.3)  smooths  out  much  of  this  noise,  giving  a  velocity  record 
which  compares  favorably  with  the  test  data.  It  should  be  noted  that  the  test 


53 


data  velocity  record  does  not  return  to  zero  as  is  physically  required.  Thtis, 
only  the  initial  velocity  records  can  be  compared.  In  the  integration  of  the 
velocity  data  to  give  the  deflection  record  (Figtire  4.4),  all  evidence  of 
high-frequency  oscillation  has  been  smoothed  out.  However,  it  can  also  be 
seen  in  this  figure  that  the  test  data  deflection  does  not  reach  a  constant 
value  before  the  end  of  the  data  record.  Thus,  any  comparisons  between 
calculated  total  deflection  and  measured  total  deflection  are  suspect. 

This  problem  In  late -time  deflection  measurement  Is  unfortunate ,  but  It 
Illustrates  a  very  common  problem  In  the  measurement  of  In- structure  shock 
with  accelerometers.  The  range  of  an  accelerometer  must  be  set  high  enou^  to 
capture  the  high  Initial  accelerations  and  as  a  result  lacks  the  sensitivity 
to  record  late -time  accelerations.  These  late- time  accelerations  are  what 
account  for  the  overall  rigid-body  motxon  of  the  structure,  and  hence  total 
deflection. 

Given  the  problem  of  high-frequency  noise  In  the  calculated  acceleration 
record  and  lack  of  sensitivity  to  measure  late  time  accelerations  and  hence 
total  deflections  in  the  test  data,  It  Is  difficult  to  make  meaningful 
comparisons  between  the  calculated  and  measured  motion- time  histories. 

However,  the  generation  of  a  shock  spectrum  (Figure  4.5)  is  driven  by  the 
velocity  time  history,  which  has  the  fewest  problems  In  both  cases.  Calcu¬ 
lated  spectral  acceleration  and  velocity  compare  quite  well,  while  deflection 
comparisons  must  be  made  with  care  due  to  the  above-mentioned  problems  in  the 
test  data. 

Calculated  vertical  acceleration  (Figure  4.6)  at  the  mld-Qoor  does  not 
^how  much  high  frequency  noise  and  compares  favorably  with  the  test  data. 
Vertical  velocity  (Figure  4.7)  and  initial  deflection  (Figure  4.8)  show  fair 
comparisons.  The  vertical  shock  spectra  (Figure  4.9)  show  favorable  compari¬ 
sons  for  peak  accelerations  and  velocities  and  a  reasonable  value  for  deflec¬ 
tion. 
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Figure  4.1.  CONUEB  2,  STABLE  Coarse  Finite- Element  Grid. 
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10 


CONWEB  2,  hid- Floor  Horizontal  Shock  Spectra,  Coa 
STABLE  Analysis  vs.  Test  Data. 


4.3.2  CONUEB  2  Fine  Grid  Analysis 

The  problem  of  high  frequency  noise  In  the  calculated  acceleration 
record  was  thought  to  Indicate  a  problem  in  the  modeling  of  this  structure 
using  a  coarse  grid.  A  finer  grid  analysis  was  therefore  conducted. 

Figure  4.10  shows  the  finer  grid  used  In  this  analysis.  All  other  parameters 
remained  the  some  as  the  coarse  grid  analysis. 

Results  of  a  10 •msec  calculation  of  the  response  of  the  CONUEB  2 
structure  using  the  fine  grid  are  shown  in  Figure  4.11  to  Figure  4,16.  More, 
not  less,  high  frequency  noise  is  evi^nt  in  both  the  horizontal  and  vertical 
acceleration  records.  Making  the  grid  finer  did  not  correct  this  problem,  and 
another  explanation  must  be  found  for  this  phenomenon.  Again,  the  velocity 
and  deflection  records  showed  less  noise,  and  in  fact,  showed  a  good  correla¬ 
tion  to  the  test  data. 

In  the  calculation  of  this  fine  grid  analysis,  another  problem  with 
STABLE  was  observed.  This  relatively  simple  problem  required  16  hours  of  run 
time  on  a  20-MHz,  386  personal  computer.  The  ultimate  aim  of  this  research 
project  is  the  qualification  of  an  in-structure  shock  analysis  tool,  which 
will  operate  on  Just  such  a  machine.  Run  times  such  as  this  are  orders  of 
oiagnltude  too  large  for  this  application. 
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Figure  4.10.  CONUEB  2,  STABLE  Fine 
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Figure  4.11.  CONVEB  2,  Mld*Floor  Horizontel  Acceleration,  Fine  Grid 
STABLE  Analysis  vs.  Test  Data. 
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Figure  4.12.  CO.SVEB  2.  Mid-Floor  Horizontal  Velocity,  Fine  Grid  STABLE 
Analysis  vs.  Test  Data. 
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Figure  4.13.  CONWEB  2.  Mid-Floor  Horizontal  Deflection,  Fine  Grid  STABLE 
Analysis  vs.  Test  Data. 
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Figure  4.14.  CONWEB  2.  Mid-Floor  Vertical  Acceleration.  Fine  Grid  STABLE 
Analysis  vs.  Test  Data. 


4.3.3  Soil -Spring  Fine  Grid  Analysis 

In  a  final  attempt  to  determine  the  cause  of  excessive  noise  In  the' 
calculated  acceleration  records,  the  fine  grid  STABLE  calculation  of  CONWEB  2 
was  repeated.  In  this  analysis,  the  SMI  at  the  back  face  of  the  structure  was 
replaced  by  equivalent , soil  springs.  It  was  hypothesized  that  problems  in  the 
calculation  of  the  relative  motion  of  the  structure  and  the  soil  on  a  non- 
loaded  interface  could  cause  fictitious  loads  to  be  generated.  These  false 
loads  could  account  for  high  frequency  noise  in  the  acceleration  record. 

The  soil -spring  properties  included  in  this  calculation  were  generated 
by  assuming  an  elastic  soil  response.  A  typical  soil  modulus  of  elasticity  of 
1,000  psi  was  assumed  based  on  recommendations  by  Bowles  [12].  Soil-spring 
stiffness  was  calculateo  based  on  an  effective  depth  of  elastic  response  of  18 
inches  at  the  back  of  the  structure.  This  is  a  very  crude  model  of  the 
complex  SMI  at  the  back  of  the  structure,  but  should  be  sufficient  to  examine 
the  problem  of  SMI -induced  instability. 

Results  of  a  20 -msec  calculation  of  the  response  of  the  CONWEB  2 
structure  using  the  fine  grid  are  shown  in  Figure  4.17  through  Figure  4.22. 
High  frequency  noise  is  evident  in  both  the  horizontal  and  vertic  accelera¬ 
tion  records.  There  was  less  noise  than  in  the  fine  grid  analysis,  but  more 
noise  chan  in  the  coarse  grid  analysis.  Removing  Che  SMI  from  Che  model  did 
not  correct  this  problem.  This  implies  Chat  the  phenomenon  is  related  to  some 
other  problem  in  the  analysis.  Again,  Che  velocity  and  deflection  records 
showed  less  noise  and,  in  fact,  showed  a  very  good  correlation  to  the  test 
data. 
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Figure  4.17.  CONWEB  2,  Mid-Floor  Horizontal  Acceleration,  Soil-Spring 
Fine  Grid  STABLE  Analysis  vs.  Test  Data. 
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Figure  4.18.  CONWEB  2,  Mid-Floor  Horizontal  Velocity,  Soil-Spring  Fine 
Grid  STABLE  Analysis  vs.  Test  Data. 
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Figure  4.19.  CONWEB  2,  Mid- Floor  Hori 
Grid  STABLE  Analysis  vs. 
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Figure  4.20.  CONWEB  2.  Mid-Floor  Verti 
Grid  STABLE  Analysis  vs. 
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4.4  STABLE  Analysis  Sunnnary 

The  evaluation  of  STABLE  as  an  in- structure  shock  analysis  technique 
leads  to  the  rejection  of  this  program  from  further  consideration.  This 
rejection  is  not  based  on  the  problem  of  high  frequency  noise  in  the  accelera¬ 
tion  records  that  plagued  this  analysis.  Favorable  comparisons  of  the 
calculations  results  to  CONUEB  2  test  data  indicate  that  the  progra^a  ia 
capable  of  handling  this  type  of  application.  Modification  of  the  program  to 
eliminate  the  proble.m  of  high  frequency  noise  in  the  acceleration  records  is 
quite  possible.  The  basis  for  the  rejection  of  the  application  of  STABLE  to 
in- structure  analysis  is  excessive  run  times  on  the  target  computer.  The  fine 
grid  analysis  required  api^roxlmately  16  hpurs  to  complete  on  a  20  MHz,  386 
personal  computer. 

Excessive  run  times  for  STABLE  are  directly  related  to  its  formulation 
as  an  implicit  finite-element  analysis  program.  As  discussed,  above  the 
implicit  formulation  refers  to  the  solution  method  used  to  solve  the  equations 
of  motion  as  the  problem  proceeds  through  time.  Using  a  large  time  step  is 
the  main  advantage  of  such  a  formulation;  however,  with  each  time  step,  the 
solution  of  the  equations  of  motion  requires  a  relatively  large  amount  6f 
calculatlonal  effort  and,  hence,  computer  time.  Unfortunately,  the  applica¬ 
tion  of  STABLE  to  in- structure  shock  analysis  requires  a  small  time  step  and  a 
large  number  of  iterations  to  capture  the  high  frequency  response  of  the 
structure  and  to  keep  the  SMI  calculation  from  going  xinstable.  As  a  result, 
you  have  a  combination  of  the  worst  of  both  worlds;  that  is,  the  requirement 
for  many  short  time  steps  that  require  a  large  amount  of  computer  time  for 
each  calculation.  Modifications  are  possible  that  could  increase  the  computa¬ 
tional  efficiency  of  STABLE,  but  it  is  felt  that  such  modifications  could  not 
overcome  such  a  fundamental  limitation. 
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CHAPTER  5 

EXPUCIT  FINITE-ELEMENT  TECHNIQUE 


5.1  ISSV3  Program 

With  Che  rejection  of  STABLE,  a  new  candidate  program  was  selected  for 
development  as  an  in-structure  shock  analysis  tool.  This  program,  ISSV3,  is 
an  explicit  type  finite-element  analysis  program  in  which  the  time  step  must 
be  small  to  assure  stability,  but  the  computational  effort  Is  smaller  at  each 
step  or  Iteration.  As  will  be  seen  below,  the  result  Is  a  program  with 
orders-of-magnltude  shorter  run  times. 

ISSV3  was  developed  by  Robert  E.  Walker,  James  L.  Drake,  William  L. 

Boyt,  and  Thomas  R.  Slawson  of  Applied  Research  Associates,  Inc.,  Southern 
Division.  The  work  was  conducted  under  WES  contract  DACA39-90-0041  for  Dr.  J. 
P.  Balsara  and  under  supervision  of  Che  author  of  this  report,  Mr.  Richard 
Dove.  Included  in  this  work  was  a  report  written  by  Walker,  Drake,  Boyt,  and 
Slawson  (5].  The  following  brief  discussion  of  the  formulation  of  ISSV3  draws 
extensively  from  information  in  this  report  but  will  foctis  on  the  application 
of  the  program  to  in- structure  shock  analysis  of  underground  structures. 

After  this  discussion,  ISSV3  will  be  applied  to  In-structure  shock  analysis  of 
the  CONVEB  test  series,  and  the  results  will  be  >^valuated  against  test  data 
and  compared  to  ocher  calculation  techniques. 

5.2  ISSV3  Program  Formulation 

ISSV3  is  a  two-dimensional,  lumped-parameter  beam  model,  explicit 
finite-element  analysis  progrem.  It  includes  nonlinear  structural  behavior, 
SMI ,  and  the  calculation  of  the  ground  shock  or  airblast  loads .  Given  the 
distance  of  a  bomb  from  a  structure  and  backfill  properties,  the  program 
calculates  the  free- field  environment.  These  free- field  loads  are.  converted 
to  structural  nodal  loads  by  an  SMI  model. 

The  central  finite-difference  direct  integration  method  is  used  to  solve 
the  equations  of  motion  as  the  solution  progresses  through  time.  As  such,  the 
displacements  and  rotations  at  a  given  time  step  are  calculated  from  the 
displacements  and  rotations  from  the  prior  time  step.  From  these  calculations 
the  nodal  displacements  and  rotations  are  \ised  to  calculate  element  end 
moments,  shears,  and  thrusts.  These  forces  are  lued  to  update  the  equations 
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of  motion  for  the  next  calculatlonal  step.  Damping  is  included  as  an  internal 
force  for  each  element  not  as  a  global  term.  One  of'  the  main  advantages  of 
this  solution  scheme  is  that  the  stiffness  and  mass  matrices  do  not  need  to  be 
reassembled  for  each  time  step,  greatly  decreasing  calculatlonal  effort.  This 
solution  method  is  widely  covered  in  the  literature  and  those  Interested 
should  refer  to  [16]  by  Bathe  and  Wilson. 

Program  output  includes  nodal  displacements,  velocities,  and  accelera¬ 
tions.  This  motion  is  the  in-structure  shock  response  for  each  node  in  the 
structure.  A  graphical  representation  of  the  rigid-body  motion  and  the  shape 
of  the  deformed  grid  is  available.  In  addition  to  the  in- structure  shock 
information  generated,  structural  response  data  such  as  element  thxrusts, 
shears,  moments,  and  strains  are  also  output. 

5.2.1  Free-Field  Load  Generation 

the  equations  used  in  ISSV3  to  characterize  the  free-field  environment 
are  essentially  those  found  in  the  TM5  855-1  [2).  As  such,  the  stress  and 
particle  velocities  at  a  given  point  in  the  free  field  are  described  by  an 
exponential  decay  time  history. 

Equations  (18)  and  (19)  are  used  in  ISSV3  for  the  pressure- time  and 
velocity- time  histories  at  a  given  point.  As  the  range  from 

P(t)  «  tzO 

P(t)  •  Free-field  stress.  (18) 

P,  ■  PeeUc  free-field  stress, 
t,  ■  Time  of  arrival, 
a  ■  Time  constant  (1.0)  . 

the  bomb  increases,  the  amplitude  of  each  history  decreases  and  its  duration 
increases. 

V{t)  »  V,(l-pt/t,)e’'*'V  tzO 

v{t)  "  Free-field  velocity.  (19) 

Vg  ■  Peak  radial  partical  velocity. 

t,  ■  Time  of  arrival. 

P  •  Time  constant,  (a/8.5)  . 

It  can  be  seen  in  these  equations  Chat  the  amplitude  and  duration  of  both 
pressure- time  and  velocity- time  histories  are  dependent  on  the  time  of  arrival 
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(t,).  In  the  default  mode  ISSV3  calculates  a  t,  is  dependent  only  on  range 
and  average  seismic  velocity.  It  has  been  found,  Hayes  [1] ,  that  this  may  be 
an  oversimplification.  Equation  (20)  shows  that  ISSV3  can  also  use  a  non¬ 
linear  arrival  time  which  Is  a  function  of  peak  stress,  peak  free -field 
velocity  and  backfill  properties. 

M 

C.  -  f  Ci  -  C+SV, 

•i 

e„  ■  Time  of  arrival  (sec) .  /20) 

R  ■  Range  (ft) . 

Cl  ■  Loading  Have  Velocity  (ft/sec)  . 
c  a  Seismic  velocity  (ft/aec) . 

S  ■  Backfill  Variable  «  l/(Air  Filled  Voids) 

V,  ■  Peak  Free  Field  Velocity  (ft/aec)  . 

Time  constants  a  and  (I  are  generally  taken  to  be  1.0  and  1/8.5,  respec¬ 
tively.  Po  and  Vq  are  the  peak  presstire  and  particle  velocity  at  the  point  of 
Interest  and  are  calculated  using  (21)  and  (22).  Pq  is  dependent  on  the 
backfill  properties  of  mass  density,  seismic  velocity,  and  attenuation  coeffi¬ 
cient,  as  well  as  the  rdnge,  bomb  weight,  and  coupling  factor. 


Po  •  160fpc(.j^)- 


P,  ■  Peak  free- field  stress  (psi)  . 
f  ■  Coupling  factor, 
n  ■  Attenuation  coefficient. 
p  ■  Mass  density  (Ib-aec* / ft*)  . 
c  ■  Seismic  velocity  (ft/aec) 

R  •  Range  (ft) . 

W  ■  Charge  weight  (lb)  . 


(21) 


Vg  Is  dependent  only  on  the  attenuation  coefficient  and  the  range,  bomb 
weight,  and  coupling  factor. 


Vg  ■  Peak  particle  velocity  (ft/aec) . 
f  ■  Coupling  factor, 
ne  Attenuation  coefficient. 

R  m  Range  (ft) . 

tf  ■  Charge  weight  (lb)  . 


(22) 


These  parameters  were  discussed  In  some  detail  in  the  TMS  855-1  In-structure 
shock  analysis  aiection  of  this  report.  As  mentioned  above,  the  characteriza¬ 
tion  of  the  free- field  environment  using  these  equations  is  thought  to  be  an 
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overslmpliftcatlon.  There  Is  research  underway  to  address  this  problem,  and 
when  better  models  are  available,  they  will  be  incorporated  into  ISSV3. 

5.2.2  Structure  Media  Interaction 

Free -field  loads  muist  be  applied  to  the  structure  through  SHI.  The  SMI 
models  used  in  ISSV3  are  similar  to  the  ones  used  in  the  STABLE  program 
discussed  above.  Interface  pressures  are  dependent  on  the  relative  motion  of 
the  backfill  and  the  structure  as  well  as  the  backfill  properties.  Equations 
(23)  and  (24)  are  used  to  calculate  the  interface ■ normal  and  shear  stresses  at 
each  node  on  the  outside  of  the  structure. 


®i  ■  ♦  PCp{V^-i^)  i  0 

o^  ■  SMI  interface  atreaa. 

a„  ■  Nortoal  ground  shock  stress. 

p  >  Soil  mass  density. 

Cp  ■  Free-field  pompression  wave  speed. 

Fp  ■  Normal  free-field  soil  velocity  vector. 
X„  •  Normal  structure  velocity  vector. 


ti  -  pc,lV^-it)  s  or  Y 

■  SMI  interface  shear  stress, 
a„  ■  Normal  ground  shock  stress, 
p  •  Soil  mass  density. 

c,  ■  Free-field  shear  wavespeed.  (f^) 

Vj  ■  Tangential  free-field  soil  velocity  vector. 

Xf  ■  Tangential  structure  velocity  vector. 

W  ■  Gravity  stress, 
y  ■  Ultimate  shear  stress. 
p  ■  Coefficent  of  Friction. 


It  should  be  noted  that  the  free-field  and  SMI  models  used  in  1S5V3  do  not 
take  into  account  the  fact  that  the  structure  perturbs  the  free-field  stress. 
In  other  words,,  the  free-field  stress  is  calculated  as  if  the  structure  does 
not  exist,  and  is  then  applied  to  the  structure  through  SMI.  This  is  thought 
CO  be  a  conservative  assumption,  and  might  best  be  considered  on  a  case-by- 
ase  basis. 

5.2.3  Internal  Forces 

A  standard  beaip  element  is  Che  element  type  used  in  ISSV3.  The  internal 
element  forces  are  calculated  using  Equations  (25),  (26),  and  (27).  The  new 


thmsts  at  each  tine  step  are  calculated  from  the  old  thrusts  plus  the 
difference  due  to  relative  displacement  of  the  end  nodes  during  the  time  step. 

I’’— •  ^ 

T  ■  Element  Thrust. 

S  ■  Element  Youngs  Modulus. 

A  ■  Element  Area. 

L  ■  Element  Length. 

The  new  moments  are  calculated  from  the  old  moments  plus  the  moments  due  to 
the  end  rotations. 


M  ■  Element  End  Moments. 

B  •  Youngs  Modulus. 

I  ■  Moment  of  Inertia. 

L  ■  Element  Length. 

■  Mear  ana  Far  End  Element  End  Rotation. 


(26) 


The  new  shears  are  developed  from  the  unbalanced  new  moments. 

Vi  -  -V, 

(27) 

V  •  Element  End  Shears. 

M  ■  Element  Moments. 

L  •  Element  Length.  ■ 

When  the  ultimate  values  of  the  thrusts  and  moments  are  reached  plastic 
deformation  occurs.  These  ultimate  values  are  calculated  prior  to  the 
analysis  based  on  Che  section  properties.  A  moment  thrust  diagram  is  con* 
structed,  and  a  balanced  pair  of  ultimate  values  are  selected. 

'  Internal  damping  is  included  as  an  internal  force  to  reduce  high 
frequency  oscillations.  IT>e  interaction  of  the  structure  and  the  surrounding 
soil  lead  to  situations  were  the  structure  moves  away  from  the  soil  or 
cavitation  cakes  place.  In  such  a  situation  the  damping  usually  provided  by 
the  soil  ceases  to  exist.  Internal  damping  forces  are  calculated  via  Equa* 
tlons  (28)  and  (29)  to  damp  out  numerical  oscillations  associated  with  the 
axial  and  double  curvature  flexural  mode  of  each  element.  These  forces  are 
transformed  from  Che  local  to  global  coordinates  and  summed  into  Che  internal 
force  vector. 
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CX  •  D  C.  Cg  - 

at 

CX  ■  Axial  Dam>i-n9  Fozca. 

D  >  Input  Damping  Ratio. 

Cg  •  Critical  Damping  Ratio. 
AL  ■  Change  in  Length. 

At  •  Time  Step. 

A  ■  Element  Area. 

E  ■  Youngfs  Modules. 
p  >  Density. 


O#  »  Z>  C,{ 


At  ‘ 


C, 


4 


N 


ML  i2LL^ 

L  24 


CM  •  Dancing  Force. 

D  ■  Input  Damping  Ratio. 

Cg  ■  Critical  Damping  Ratio. 
9  ■  Element  End  Rotation. 

At  »  Time  Step. 

E  ■  Youngs  Modulus. 

I  •Moment  of  Inertia. 

L  •  Element  Length. 

A  a  Element  Area. 
p  a  Density. 


A  fiLI) 


(29) 


3,2.4  Equilibriun  Equations 

The  equilibriun  equations  solved  at  each  tine  step  are  cpnstructed  fron  the 
internal  and  external  forces  transfomed  into  the  global  coordinate  systea  as 
shorn  in  Equation  (30).  As  discussed  before,  the  central  difference  integra¬ 
tion  method  is  used  to  solve  this  equation  for  notions  as  the  solution 
progresses  through  time. 

X  -  {Fg-Fj)/M 

Jt  a  Nodal  Acceleration.  (30) 

Fg  •  External  Forces. 

Fj  •  Internal  Forces. 

M  a  Nodal  mass. 


5.3  ISSV3  conwe: 

Oevelopneq 
requires  some  pi[ 
selected  from  a 
path  is  selected 
condition,  Usin 


B  Analyses 

t  of  section  properties  for  inputting  into  the  ISSV3  prograa 
eanalysis.  The  ultimate  moment  and  thrust  capacities  must  be 
calculated  moment- thrust  diagraa  for  each  section.  A  load 
giving  a  moment- thrust  velue  Just  below  the  balanced 
g  the  material  properties  for  the  CONVEB  test  series,  section 
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prcpercles  w«re  generated  for  each  conponent.  These  section  properties  are 
shown  in  Table  S . 1 . 


Table  S.l.  ISSV3,  CONVEB  Section  Properties. 


Test 

Coaponent 

E 

pai 

P 

Ib/ft^ 

I 

in* 

A 

in* 

Mult 
in- lb 

Tult 

lb 

CONVEB  1,  3,  4 
Test  Vail 

4E6 

150 

B 

1.42E4 

4.58E3 

CONVEB  2 

Test  Vail 

4E6 

150 

32.0 

8.6 

S.90E4 

1.30E4 

Reaction 

Structure 

4E6 

150 

83 

11 

1.8E5 

2.4E4 

5.3.1  CONVEB  1  Analysis 

The  finite 'eleiien'.:  grid  for  the  ISSV3  ln>struct\ure  analysis  of  CONVEB  1 
is  shorn  in  Figure  S.l.  This  saae  grid  was  u4ed  for  the  calculation  of  CONVEB 
3  and  4,  with  nodificstions  to  the  backfill  properties.  Section  properties 
used  were  shown  above  in  Table  S.l,  and  backfill  properties  were  shown  in 
Table  2.1. 

Results  froh  the  CONVEB  1  analysis  are  shown  in  Figure  S.2  to 
Figure  S.23.  It  is  iaoediately  obvious  on  exaaination  of  the  interface 
pressure  coaparison  that  the  calculated  time  of  arrival  (t,)  of  the  load  at 
the  face  of  the  structure  lags  the  actual  value  by  about  4  asec.  This 
phenoaenon  was  observed  throughout  the  CONVEB  calculations.  Interestingly 
enough,  the  shape  and  aagnitude  of  the  calculated  interface  stress  coapare 
very  favorably  with  the  test  data,  in  spits  of  the  strong  relationship  between 
t,  and  free* field  stress. 

As  a  result  of  the  tiae  lag  in  the  loading,  all  the  data  eo8n>«rlsona  for 
structural  response  also  lag.  Coaparisons  tc  test  data  of  calculated  acceler* 
ation,  velocity,  deflection  histories,  and  shock  spectfr  are  shown  in 
Figure  S.4  to  Figure  S.23.  Horizontal  acceleration,  velocity,  and  deflection 
histories  were  reproduced  reasonably  well  by  the  ISSV3  analysis.  Exaaination 
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of  the  late-tlae  horizontal  velocity  records  shows  that  ISSV3  underpredicts 
the  late-tlae  velocity.  This  is  thought  to  stea  froa  the  inability  of  the 
current  free-field  aodel  to  capture  the  late-tiae  flow  of  a  clay-type  aateri- 
al,  a  problea  which  should  be  solved  In  later  versions  of  the  program.  Also, 
as  discussed  in  the  STABLE  analysis  section,  there  Is  great  difficulty 
aeasuring  rigid-body  displacements  with  acceleroaeters,  due  to  the  low  level 
of  acceleration  for  these  late-tiae  actions. 

Vertical  acceleration,  velocity,  and  deflection  histories  were  also 
reproduced  reasonably  well  by  the  ISSV3  analysis.  Due  to  the  snich  saaller 
aagni tildes  of  the  vertical  velocities  and  deflections,  uncertainties  in 
aeasureaent  of  late-tiae  aotions  have  a  larger  relative  effect  than  in  the 
horizontal  measurements.  Given  the  uncertainties  of  these  aeasureaents ,  the 
magnitudes  calculated  for  vertical  velocity  and  deflection  by  ISSV3  are 
reasonable. 

Shock  spectra  were  generated  from  the  velocity  histories  for  each 
position  in  Che  structure.  These  spectra  show  a  good  correlation  between  the 
calculated  values  and  spectra  generated  for  test  data.  The  only  consistent 
underprediction  was  for  Che  rigid-body  deflection  as  seen  on  the  low  frequency 
portion  of  the  spectra.  This  underpredlction  is  thought  to  be  due  Co  a 
coabinaclon  of  the  problea  with  the  aodellng  of  clay  backfill  in  the  analysis 
and  the  problea  of  aeasuring  late-tiae  aotions  with  acceleroaeters  . 
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Figure  5.11.  CONWEB  1,  Mid- Floor  Horizontal  Shock  Spectra,  ISSV3  Node  9 
Analysis  vs.  Test  Data. 
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Figure  5.18. 
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Figure  5vl9. 
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Figure  5.22.  CONWEB  1.  Floor  Vertical  Deflection,  ISSV3  Node  8, 
Analysis  vs.  Test  Data. 
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Figure  3.23.  CONUEB  1,  Floor  Vertical  Shock  Spectra,  ISSV3  Node  8, 
Analyi is  vs.  Test  Data. 
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5.3.2  CONWEB  2  Analysis 

As  discussed  above.  CONWEB  2  was  a  test  with  the  same  clay  backfill  as 
in  CONWEB  1.  The  test  wall  was  twice  as  thick  as  in  the  other  tests,  with  an 
L/d  ratio  of  5.  The  finite-element  grid  for  the  ISSV3  in-structure  shock 
analysis  of  CONWEB  2  is  shown  in  Figure  5.24.  This  grid  is  different  f ro.  .  the 
one  used  for  the  calculation  of  CONWEB  1,  3,  and  4,  due  to  the  different  test 
wall.  Section  properties  used  in  this  an.ily3is  were  shorn  in  previously 
Table  5.1. 

Results  from  the  CONWEB  2  analysis  are  shown  in  Figure  5.25  to 
Figure  5.46.  As  in  CONWEB  1.  examination  of  the  interface  pressures  show  that 
the  calculated  time  of  arrival  (t,)  of  the  load  at  the  face  of  the  structure 
lags  the  actual  value  by  about  4  msec.  Again,  the  shape  and  magnitude  of  the 
calculated  interface  stress  compare  very  favorably  with  the  test  data,  in 
spite  of  the  strong  relationship  between  t^^  and  free-field  stress. 

As  before,  all  structural  response  data  also  lag  by  about  4  msec. 
Comparisons  to  test  data  of  calculated  acceleration,  velocity,  deflection 
histories,  and  shock  spectra  are  shown  in  Figure  5.27  to  Figure  5.46. 
Horizontal  acceleration,  velocity,  and  deflection  histories  were  reproduced 
reasonably  well  by  Che  ISSV3  analysis.  The  front  wall  response  (Figure  5.27 
to  Figure  5.30)  shows  an  overpredictlon  of  velocity  and  deflection  response. 
Even  here,  examination  of  the  late  time  horizontal  velocity  records  show  that 
ISSV3  under  predicts  the  late-cime  velocity  due  Co  late-time  clay  backfill 
effects.  On  Che  whole,  however,  the  calculated  horizontal  in-structure  shock 
response  compares  very  favorably  with  Che  test  data  throughout  the  structure, 
given  Che  difficulty  of  measuring  rigid-body  displacements  with  acceler¬ 
ometers.  .  CONWEB  2  vertical  acceleration,  velocity,  and  deflection  histories 
were  also  reproduced  reasonably  well  by  the  ISSV3  analysis.  As  in  CONWEB  1 
Che  small  magnitudes  of  the  vertical  velocities  and  deflections  cause  C^e 
uncertainties  in  measurement  of  late- time  motions  to  have  a  large  relative 
effect.  Given  Che  uncertainties  of  these  measurenents ,  the  magnitudes  calcu¬ 
lated  for  vertical  velocity  and  deflection  by  ISSV3  are  reasonable. 

Shock  spectra  were  generated  from  Che  velocity  histories  for  each 
position  in  the  structure.  These  spectra  show  a  very  good  correlation  between 
Che  calculated  values  and  spectra  generated  for  test  data.  As  in  CONWEB  1  Che. 
only  consistent  underprediction  was  for  the  rigid-body  deflection  as  seen  on 
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Che  low-frequency  portion  of  the  spectra.  Again,  this  underprediption  is 
thought  to  be  due  to  a  combination  of  the  problem  with  the  modeling  of  clay 
backfill  in  the  analysis,  and  the  problem  of  measuring  late-time  motions  with 
accelerometers.  The  front  wall  response  was  slightly  over-predicted,  which 
implies  that  the  model  of  the  wall  was  softer  than  in  the  acCtial  test.  The 
modeling  of  the  front  wall  as  a  one-way  slab,  neglecting  all  two-way  effects, 
could  account  for  this  relatively  minor  effect.  Also,  the  fact  that  this 
phenomenon  did  not  occur  in  the  thinner  test  wall  in  CONWEB  1  seems  to  imply 
that  the  difference  is  due  to  the  relative  importance  of  shear  response  for 
the  thicker  CONUEB  2  wall.  However,  it  will  be  seen  later  that  the  response 
of  the  thin  walls  in  CONWEB  3  and  4  was  overpredicted  -iiich  does  not  support 
this  hypothesis.  Shear  response  is  not  modeled  in  v  current  version  of  the 
ISSV3  program. 
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.24.  CONVEB  2,  ISSV3  Finite -ElenenC  Grid  Superiapoeed  on 
Structure  Cross  Section. 
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CONUEB  2,  Mid-Floor  Horizontal  Shock  Spectra,  ISSV3  Node  9 
Analysis  vs.  Test  Data. 
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Figure  5.42.  CONVEB  2,  Floor  Horizontal  Shock  Spectra,  ISSV3  Node  8, 
Analysis  vs.  Test  Data. 
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Figur«  5.43.  CONWEB  2,  Floor  Vertical  Acceleration,  ISSV3  Node  8, 
Analysis  vs.  Test  Data. 
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Figure  5.44.  CONWEB  2,  Floor  Vertical  Velocity,  ISSV3  Node  8, 
Analysis  vs.  Test  Data. 
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.45.  CONWEB  2,  Floor  Vertical  Deflection,  1SSV3  Node  8, 
Analysis  vs.  Test  Data. 
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.46.  CONWEB  2,  Floor  Vertical  Shock  Spectra,  1SSV3  Node  8 
Analysis  vs.  Test  Data. 
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5.3.3  COWEB  3  Analysis 

CONWEB  3  was  a  test  with  the  sane  structure  as  in  CONWEB  1  and  4,  the 
test  wall  had  the  sane  thickness,  with  an  L/d  ratio  of  10.  In  CONWEB  3  a  sand 
backfill  was  used.  The  finite-elenent  grid  for  the  ISSV3  in-structure 
analysis  of  CONWEB  3  is  shown  in  Figure  5.1.  This  savie  grid  was  used  for  the 
calculation  of  CONWEB  1  and  4,  with  modifications  to  the  backfill  properties. 
Section  properties  used  were  shown  previously  in  Table  5.1,  and  backfill 
properties  were  shoyn  in  Table  2.1. 

CONWEB  3  analysis  results  are  presented  in  Figure  5.47  to  Figure  5.68. 

As  in  CCNWEB  1,  comparison  of  the  interface  pressures  once  again  shows  that 
the  calculated  time  of  arrival  (t,)  of  the  load  at  the  face  of  the  structure 
lags  the  actual  value  by  about  4  msec.  Again,  the  shape  and  magnitude  of  the 
calculated  Interface  stress  at  the  center  of  the  wall  (Figure  5.48)  compare 
favorably  with  the  test  data.  Figure  5.47  shows  the  interface  pressure 
measured  on  the  lower  part  of  the  test  wall  and  illustrates  an  interesting 
effect  of  the  clynamic  loading  of  buried  structures  in  sands.  This  effect, 
known  as  soil  arching,  results  from  the  redistribution  of  interface  loads  from 
the  center  of  a  loaded  span  to  the  supports.  The  exirrent  free -field  model  in 
ISSV3  does  not  include  this  effect,  but  this  should  have  only  minor  impact  on 
the  overall  In-structure  shock  response. 

As  before,  all  the  data  comparisons  for  structural  response  also  lag  by 
about  4  msec  due  to  the  loading  t, .  Comparisons  to  test  data  of  calculated 
acceleration,  velocity,  deflection  histories,  and  shock  spectra  are  shown  in 
Figure  5.49  to  Figure  5.68.  Horizontal  acceleration,  velocity,  and  deflection 
histories  were  reproduced  reasonably  well  by  the  ISSV3  analysis.  There  is  no 
evidence  of  the  underprediction  of  late  time  velocity  in  this  sand ' backfill 
test.  This  lends  credence  to  the  hypothesis  that  the  late-time  velocity 
underpredictions  in  CONWEB  1  and  2  were  due  to  flow  effects  in  the  clay 
backfill. 

AS  in  CONWEB  2,  the  front  wall  response  (Figure  5.49  to  Figure  5.52) 
shows  un  overprediction  of  velocity  and  deflection  response.  Given  that  the 
accelerrtion  history  compares  very  well  at  this  location,  the  overprediction 
of  velocity  and  deflection  implies  that  the  model  of  the  front  wall  is  softer 
than  the  acttial  wall.  On  the  whole,  however,,  the  calculated  horizontal  in- 


structure  shock  response  compares  favorably  with  test  data  throughout  the 
structure . 

CONWEB  3  vertical  acceleration,  velocity  and  deflection  histories  were 
also  reproduced  reasonably  well  by  the  ISSV3  analysis.  The  magnitudes  of  the 
vertici.1  velocities  and  deflections  were  even  smaller  than  in  CONWEB  1  and  2, 
causing  the  uncertainties  in  measurement  of  late-time  motions  to  have  an  even 
larger  relative  effect.  Given  the  uncertainties  of  these  measurements,  the 
magnitudes  calculated  for  vertical  velocity  and  deflection  by  ISSV3  are 
reasonable. 

Shock  spectra  were  generated  from  the  velocity  histories  for  each 
position  in  the  structure.  These  spectra  show  a  good  correlation  between  the 
calculated  values  and  spectra  ganerated  for  test  data.  Rigid-body  deflectioti 
as  seen  on  the  low-frequency  portion  of  the  spectra  was  not  uhderprcdicted  for 
this  sand  backfill  test  as  it  was  fox*  clay  backfill  tests  in  CONWEB  1  and  2. 
This  indicates  that  the  underprediction  in  CONWEB  1  and  2  was  due  to  the 
problem  with  the  modeling  late  time- flow  effects  of  the  clay  backfill. 

As  discussed  above,  the  front  wall  response  was  slightly  overpredicted, 
a  result  that  implies  that  the  wall  model  was  softer  than  the  actual  wall. 

The  modeling  of  the  front  wall  as  a  one-way  slab  neglecting  all  two-way 
effects  could  account  for  this  relatively  minor  effect.  The  fact  that  this 
phenomenon  occurred  with  this  thin  test  wall  implies  that  the  difference  may 
not  be  due  to  the  relative  importance  of  shear  response  for  the  thick  vs.  thin 
wall  as  was  hypothesized  in  the  CONWEB  2  discussion. 
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Figure  5.63.  CONWEB  3,  Floor  Horizontal  Deflection,  ISSV3  Node  8, 
Analysis  vs.  Test  Data. 
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Figure  3.64.  CONVEB  3,  Floor  Horizontal  Shock  Spectra,  XSSV3  Node  8, 
Analysis  vs.  Test  Data. 
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Figure  5.67.  CONWEB  3.  Floor  Vertical  Deflection,  ISSV3  Node  8, 
Analysis  vs.  Test  Data. 


Snock  Se«ctr»  CoMoar^tdA* 
Con««0  ] 


HZ 


Figure  5.68.  CONWEB  3,  Floor  Vertical  Shock  Spectra,  ISSV3  Node  8, 
Analysis  vs.  Test  Data. 
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5.3.4  CONWEB  4  Analysis 

CONWEB  '  was  a  test  with  the  same  structure  as  in  CONWEB  1  and  3.  The 
backfill  in  this  test  consisted  of  the  in-situ  soil  in  the  test  area,  a  very 
high  seismic  velocity,  saturated  gravely  clay  material.  The  test  wall  was  the 
same  thickness  as  in  CONWEB  1  and  3,  with  an  L/d  ratio  of  10.  The  finite- 
element  grid  for  the  ISSV3  in- structure  analysis  of  CONWEB  4  is  shown  in 
Figure  5.1.  This  same  grid  was  used  for  the  calculation  of  CONWEB  1  and  3, 
with  modifications  to  the  backfill  properties.  Section  properties  used  were 
shown  above  in  Table  5.1;  backfill  properties  were  shown  in  Table  2.1. 

Results  from  the  CONWEB  4  analysis  are  shown  in  Flg’'Xe  5.69  to 
Figure  5.90.  As  in  CONWEB  1,  examination  of  the  interface  pressure  comparison 
once  again  shows  that  the  calculated  time  of  arrival  (t,)  of  the  load  at  the 
face  of  the  structure  lags  the  actual  value.  However,  the  lag  here  is  about  2 
msec,  versus  about  4  msec  in  CONWEB  1.  The  shape  and  magnitude  of  the  calcu¬ 
lated  interface  stresses  (Figure  5.69  and  Figure  5.70)  show  a  sharper  higher 
peak  stress  than  shown  in  the  test  data.  However,  the  calculated  interface 
stress  compares  reasonably  well  with  the  test  data.  As  would  be  expected  in  a 
clay  test,  there  is  no  evidence  of  the  soil  arching  as  was  seen  in  the  CONWEB 
3  sand  backfill  test.  Soil  arching  is  a  phenomenon  associated  with  granular 
soils. 

All  the  data  comparisons ,  for  structural  response  also  lag  by  about  2 
msec  due  to  the  lag  in  loading  t,.  Comparisons  to  test  data  of  calculated 
acceleration,  velocity,  deflection  histories,  and  shock  spectra  are  shown  in 
Figure  5.71  to  Figure  5.90.  Horizontal  acceleration,  velocity,  and  deflection 
histories  were  reproduced  reasonably  well  by  the  ISSV3  analysis.  As  was  seen 
in  CONWEB  1  and  2,  there  is  some  evidence  of  the  underprediction  of  late-time 
velocity  in  this  test  due  to  late-time  flow  effects  in  the  clay  backfill. 

As  in  CONWEB  2  and  3,  the  front  wall,  response  (Figure  5.71  to 
Figure  5.74)  shows  an  overprediction  of  velocity  and  deflection  response.  On 
the  whole,  however,  the  calculated  horizontal  in-structure  shock  response 
compares  favorably  with  test  data  throughout  the  structure. 

CONWEB  4  vertical  acceleration,  velocity,  and  deflection  histories  were 
also  reproduced  reasonably  well  by  the  ISSV3  analysis.  The  magnitudes  of  the 
vertical  velocities  and  deflections  were  small  as  in  CONWEB  1,  2,  and  3;  thus, 
the  uncertainties  in  measurement  of  late-time  motions  have  a  large  t^lacive 


effect.  Given  the  uncertainties  of  these  measurements,  the  magnitudes  calcu¬ 
lated  for  vertical  velocity  and  deflection  by  ISSV3  are  reasbnaole. 

Shock  spectra  were  generated  from  the  velocity  histories  for  each 
position  in  the  structure.  These  spectra  show  a  very  good  correlation  between 
the  calculated  values  and  spectra  generated  for  test  data  on  the  front  wall. 
Horizontal  rigid-body  motion,  as  measured  at  node  9  in  the  center  of  the  floor 
vFigure  5.78),  was  underprcdicted  when  co..:p.?red  to  test  data  spectra.  This  is 
the  same  gage  (ASHIO)  that  gave  anomalous  results  in  the  prior  analyses. 
Spectra  generated  at  node  8  (Figure  5.86),  12  inches  further  back  on  the 
floor,  show  a  very  good  correlation  with  data  gathered  at  that  point  (gage 
ASHll) .  Given  the  heavy  damage  suffered  by  the  structiire  in  this  test,  gage 
failure  at  node  9  (ASHIO)  is  quite  possible. 
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Figure  5,83. 
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Figure  5.84. 
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87.  CONWEB  4,  Floor  Vertical  Acceleration.  ISSV3  Node  8 
Analysis  vs.  Test  Data. 
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Figure  5.88.  CONUEB  4,  Floor  Vertical  Velocity.  ISSV3  Node  8. 
Analysis  vs.  Test  OatA. 
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Figure  5.89.  CONUEB  4,  Floor  Vertical  Deflection,  ISSV3  Node  8, 
Analysis  vs.  Test  Data. 
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Figure  5.90.  CONWEB  4,  Floor  Vertical  Shock  Spectra,  ISSV3  Node  8, 
Analysis  vs.  Test  Data. 


129 


5.4  ISSV3  Analysis  Sujnmary 

The  Initial  evaluation  of  ISSV3  as  an  in-structure  shock  analysis 
technique  is  favorable.  Calculation  of  all  four  CONWEB  tests  with  very 
dissimilar  backfill  materials  gave  reasonable  results.  The  acceleration 
records  generated  had  forms  and  magnitudes  quite  similar  to  test  data. 

Compf  '.sons  of  velocity  and  deflection  records  were  sometrtiat  hampered  due  to 
the  difficulty  of  measuring  low-level,  late- time  accelerations  in  dynamic 
tests.  However,  available  data  compared  favorably  with  the  ISSV3- generated 
velocity  and  deflection  histories.  Shock  spectra  generated  from  the  ISSV3 
velocity  histories  also  compared  very  well  with  spectra  derived  from  test 
data. 

It  is  acknowledged  that  the  current  free-field  model  in  ISSV3  needs  to 
be  improved.  Late -time  velocities  are  sometimes  underpredicted  due  to  flow 
effects  in  clay  backfill.  Also,  ether  research  by  Hayes  [1]  and  [10]  has 
Indicated  chat  the  characterization  of  backfill  by  seismic  velocity,  density, 
and  attenuation  coefficient  may  be  overly  simplistic.  In  spite  of  these 
problems  in  the  free-field  models,  the  interface  loads  generated  by  ISSV3  for 
Che  CONUEB  tests  were  reasonable. 

The  single  greatest  strength  of  the  1SSV3  program  brought  to  light  in 
the  analysis  of  the  CONWEB  test  series  is  the  speed  of  the  calculation.  All 
of  Che  calculations  required  less  chan  four  minutes  Co  run  on  a  386,,  20-NHz 
personal  computer .  This  ia  a  marked  improvement  over  the  four  to  sixteen  hour 
run  times  for  this  application  of  the  STABLE  program.  . 
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CHAPTER  6 

COMPARISON  OF  ANALYSIS  TECHNIQUES 


With  the  success  of  the  initial  evaluation  of  the  ISSV3  program  as  an 
in-structure  shock  analysis  tool,  it  is  interesting  to  compare  the  program's 
results  with  the  other  available  tools.  Figure  6.1  to  Figure  6.13  show  direct 
comparisons  of  the  shock  spectra  generated  from  the  CONVEB  test  series  data 
and  all  of  the  means  of  analysis  investigated.  The  average  front  wall 
response  is  compared  to  the  TM5  853-1  method,  the  SDOF  wall  analysis,  the 
STABLE  results  (for  CONWEB  2),  and  the  ISSV3  results.  The  shock  response  on 
the  floor  of  each  structure  is  compared  to  the  TM5  855-1  method,  the  rigid- 
body  motion  SDOF  analysis,  the  STABLE  results  (for  CONWEB  2),  and  the  ISSV3 
results.  In  every  case,  the  results  of  the  1SSV3  analysis  are  at  least  as 
accurate  as  the  other  techniques.  With  planned  improvements  in  the  free-field 
model,  these  comparisons  are  expected  to  be  even  better  for  later  versions  of 
the  program. 

It  is  important  to  note  that  the  simplified  in-structure  shock  methods 
investigated  here  are  limited  in  a  way  that  ISSV3  is  not.  The  TM5  855-1 
method  and  the  rigid-body  SDOF  analysis  are  best  suited  to  simple,  small 
structures .  There  is  no  fundamental  reason  why  ISSV3  can  not  be  applied  to 
the  fast  and  accurate  in- structure  shock  analysi-  of  larger,  complicated, 
multifloor,  multibay  structures. 
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CONWEB  1,  Floor  Horizontal  Shock  Spectra,  ISSV3  Node  8, 
Cpmparison  of  Analysis  Techniques. 
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CONWEB  2,  Mid-Wall  Horizontal  Shock  Spectra,  isp3  Node  13, 
Comparison  of  Analysis  Techniques. 
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Figure  6.S.  CONVEB  2,  Mid- Floor  Horizontal  Shock  Spectra,  ISSV3  Node  9, 
Cooparison  of  Analysis  Techniques. 
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Figure  6.6.  CONVEB  2,  Mid-Floor  Vertical  Shock  Spectra,  ISSV3  Node  9 
Cooper ison  of  Analysis  Techniques. 


Figure  6.7.  CONVEB  2,  Floor  HorizcnCal  Shock  Spectra,  ISSV3  Node  8, 
Comparison  of  Analysis  Techniques. 
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Figure  6.8.  CONVEB  3,  Mid-Wall  Horizontal  Shock  Spectra,  ISSV3  Node  13, 
Comparison  of  Analysis  Techniques. 
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Figure  6.11.  CONWEB  4,  Mid-Wall  Horizontal  Shock  Spectra,  ISSV3  Node  13, 
CoBparison  of  Analysis  Techniques. 
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Figure  6.12. 


CONWEB  4.  Mid- Floor  Horizontal  Shock  Spectra,  ISSV3  Node  9, 
Comparison  of  Analysis  Techniques. 
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Flgura  6.13.  CONWEB  4,  Floor  Horizontal  Shock  Spectra,  1SSV3  Node  8 
Coaparlson  of  Analysis  Techniques. 


CHAPTER  7 

SUMMARY  AND  CONCLUSIONS 


7 . 1  Suamary 

The  semi -empirical  calculacional  procedure  in  TM5  855-1  was  shown  to 
give  reasonable  estimates  of  the  in-structure  shock  responses  of  the  first 
three  CONWEB  tests.  Inconsistences  in  the  gage  response  in  CONWEB  4  hampered 
comparisons  for  this  test,  but  results  appeared  to  be  reasonable.  The 
greatest  weakness  of  this  procedure  is  chat  It  was  developed  for  simple  box¬ 
like  structures.  The  assumptions  that  must  be  made  to  apply  this  procedure  to 
more  complicated  structures  have  not  been  thoroughly  evaluated.  Great  caution 
should  be  used  in  the  application  of  such  an  empirically  derived  method  to 
situations  outside  of  chose  actually  tested. 

The  results  of  the  wall  SDOF  analyses  were  relatively  inconsistent  due 
to  difficulty  in  developing  reasonable  loads  to  apply  to  Che  SDOF  model.  The 
loads  used  here  were  developed  using  Che  procedures  in  TM5  855-1  [2].  It  is 
felt  Chat  these  procedures  do  not  adequately  model  sttucture-media  interaction 
(SMI).  There  is  current  research  at  VES  to  address  these  deficiencies  in  SMI 
modeling. 

The  rigid-body  SDOF  analysis  gave  teasor<able  predictions  of  peak 
deflections  and  accelerations  for  CONVEB  1.  2,  and  3.'  Peak  velocities  were 
overpredicted  for  these  same  tests.  The  same  seemed  to  be  true  of  CONVEB  4 
when  comparisons  were  made  to  horizontal  acceleration  gage  ASHll.  However, 
CONVEB  4  predictions  in  comparison  with  gage  ASHIO  were  low  for  deflection  and 
acceleration  and  reasonable  for  velocity.  Underprediction  of  CONVEB  4  in 
comparison  with  gage  ASHIO  occurred  in  all  of  the  analysis  procedures  used, 
perhaps  lending  ciedence  to  the  possibility  of  testing  or  gage  problems.  It 
was  hypothesized  that  the  consistent  overprediction  of  peak  spectral  velocity 
in  the  rigid-body  SDOF  analysis  was  due  to  problems  in  modeling  structure 
media  interaction.  Overloading  the  structure  with  impulse  leads  to  over¬ 
prediction  of  velocity..  This  should  also  lead  to  an  overprediction  of 
deflection.  The  fact  that  deflection  predictions  are  reasonable,  indicates  a 
compensating  overrestraint  of  the  model.  The  Poisson’s  Ratios  assumed  for 
this  analysis  were  on  the  high  end  of  suggested  values  and  this  could  easily 
account  for  this  overre-straint .  Taken  together,  dll  this  illustrates  the 
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difficulty  of  analyzing  the  complex  response  of  a  burled  structure  using  a 
simple  SDOF  model.  It  should  be  noted,  however,  that  the  velocity  predictions 
were  conservative. 

The  evaluation  of  STABLE  as  an  In-structure  shock  analysis  technique  led 
to  the  rejection  of  this  program  from  further  consideration.  This  rejection 
was  not  based  on  the  problem  of  high  frequency  noise  in  the  acceleration 
records  that  plagued  this  analysis.  Favorable  comparisons  of  the  calculations 
results  to  CONVEB  2  test  data  indicate  chat  the  program  is  capable  of  handling 
this  type  of  application.,  Modification  of  the  program  Co  eliminate  Che 
problem  of  high  frequency  noise  in  the  acceleration  record  is  quite  possible 
and  could  perhaps  be  as  simple  as  Che  Inclusion  of  additional  damping  terms. 
The  basis  for  the  rejection  of  the  application  of  STABLE  Co  in*sCr\tcCuze 
analysis  is. excessive  run  times  on  the  target  computer.  The  fine  grid 
analysis  required  approximately  16  hours  to  complete  on  a  20  MHz,  386  personal 
computer.  Excessive  run  times  for  STABLE  are  directly  related  to  its  formula¬ 
tion  as  an  implicit  finite-element  analysis  program.  Modifications  are 
possible  chat  could  increase  Che  computational  efficiency  of  STABLE,  but  it  is 
felt  Chat  such  modifications  could  not  overcome  such  a  fundamental  limitation. 

The  initial  evaluation  of  ISSV3  as  an  in-structure  shock  analysis 
technique  was  favorable .  Calculation  of  all  four  CONVEB  tests  with  very  , 
different  backfill  materials  gave  very  reasonable  results.  The  acceleration 
records  generated  had  forms  and  magnitudes  quite  similar  to  test  data. 
Comparisons  of  velocity  and  deflection  records  were  somewhat  hamperedi  due  to 
Che  difficulty  of  measuring  low-level,  late-time  accelerations  in  dynamic 
tests.  However,  available  data  compared  favorably  with  the  ISSV3- generated 
velocity  and  deflection  histories.  Shock  spectra  generated  from  the  ISSV3 
velocity  histories  also  compared  very  well  with  spectra  derived  from  test 
data. 

It  is  acknowledged  chat  the  current  free-field  model  in  ISSV3  needs  to 
be  improved.  Late-time  velocities  are  sometimes  underpredicted  due  to  flow 
effects  in  clay  backfill.  Also,  ocher  research  by  Hayes  [1]  and  [10]  has 
indicated  that  Che  characterization  of  backfill  by  seismic  velocity,  density, 
and  attenuation  coefficient  may  be  overly  simplistic.  In  spite  of  these 
problems  in  the  free-field  models,  the  Interface  loads  generated  by  ISSV3  for 
the  CONVEB  tests  were  reasonable: 


The  single  greatest  strength  of  the  ISSV3  program  brou^t  to  light  in 
the  analysis  of  the  CONVEB  test  sOries  is  the  speed  of  the  calculation.  All 
of  the  calculations  required  less  than  four  minutes  to'  run  on  a  386,  20-MHz 
personal  computer.  This  compares  very  favorably  to  the  four  to  sixteen  hour 
run  times  for  this  application  of  the  STABLE  program. 

With  the  success  of  the  initial  evaluation  of  the  ISSV3  program  as  an 
in-structure  shock  analysis  tool,  a  comparison  vas  made  of  the  program's 
results  with  the  other  available  cools.  A  direct  comparison  was  made  of  the 
shock  spectra  generated  from  the  CONVEB  test  series  data  and  all  of  the  means 
of  analysis  investigated.  The  average  front  wall  response  was  compared  to  the 
TM5  855-1  method,  the  SDOF  wall  analysis,  the  STABLE  results  (for  CONVEB  2), 
and  the  ISSV3  results.  The  shock  response  on  the  floor  of  each  structure  was 
compared  to  the  TM5  855-1  method,  the  rigid-body  motion  SDOF  analysis,  the 
STABLE  results  (for  CONVEB  2),  and  the  ISSV3  results.  In  every  case,  the 
results  of  the  ISSV3  analysis  were  at  least  as  accurate  as  the  other  tech¬ 
niques.  Vith  planned  improvements  in  the  free-field  model,  these  comparisons 
are  expected  to  be  even  better  for  later  versions  of  the  program. 

It  is  important  to  note  chat  the  simplified  ln-scruct\ire  shock  methods 
investigated  here  are  limited  in  a  way  that  ISSV3  is  not.  The  TM5  855-1 
method  and  the  SDOF  analyses  are  best  limited  to  simple,  small  sCrucCiires. 
There  is  no  fundamental  reason  why  ISSV3  can  not  be  applied  to  Che  fast  and 
accurate  in-structure  shock  analysis  of  larger,  complicated,  multifloor, 
multibay  structures. 

7.2  Conclusions 

I 'SV3  has  been  shown  to  be  the  best  in-structure  shock  analysis  tool 
examined.  Extensive  comparisons  of  Che  results  of  the  in-structure  shock 
analysis  of  Che  CONVEB  test  series,  using  several  available  methods,  show -that 
1SSV3  is  fast,  accurate,  and  easy  to  use.  The  speed  with  which  ISSV3  can 
analyze  a  given  problem  will  allow  the  user  the  flexibility  to  quickly 
complete  numerous  in- structure  shock  calculations.  This  should  encourage  a 
designer  to  include  in- structure  shock  considerations  in  Che  early  design 
phase  of  %  protective  facility. 

Future  work  is  planned  to  improve  the  current  free-field  model  in 
ISSV3.  Tnq>rovemenCs  in  Che  SMI  model  are  also  being  considered  pending  Che 
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results  of  current  research  at  UES.  Further  analysis  effort  is  required  to 
validate  the  ISSV3  program  for  use  on  multifloor,  miltibay  burled  structures 
as  well  as  aboveground  structures. 
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